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ABSTRACT 

We report on our search for microlensing towards the Large Magellanic Cloud (LMC). 
Analysis of 5.7 years of photometry on 11.9 million stars in the LMC reveals 13 - 17 
microlensing events. A detailed treatment of our detection efficiency shows that this 
is significantly more than the ~ 2 to 4 events expected from lensing by known stellar 
populations. The timescales {t ) of the events range from 34 to 230 days. We estimate 
the microlensing optical depth towards the LMC from events with 2 < t < 400 days 
to be r|°° = l.2l:o.t x 10"^ with an additional 20% to 30% of systematic error. The 
spatial distribution of events is mildly inconsistent with LMC/LMC disk self- lensing, 
but is consistent with an extended lens distribution such as a Milky Way or LMC 
halo. Interpreted in the context of a Galactic dark matter halo, consisting partially 
of compact objects, a maximum likelihood analysis gives a MACHO halo fraction of 
20% for a typical halo model with a 95% confidence interval of 8% to 50%. A 100% 
MACHO halo is ruled out at the 95% C.L. for all except our most extreme halo model. 
Interpreted as a Galactic halo population, the most likely MACHO mass is between 
0.15 M0 and 0.9 Mq, depending on the halo model, and the total mass in MACHOs out 
to 50 kpc is found to be 9li'^ x 10^° M©, independent of the halo model. These results 
are marginally consistent with our previous results, but are lower by about a factor of 
two. This is mostly due to Poisson noise because with 3.4 times more exposure and 
increased sensitivity to long timescale events, we did not find the expected factor of ~ 4 
more events. Besides a larger data set, this work also includes an improved efficiency 
determination, improved likelihood analysis, and more thorough testing of systematic 
errors, especially with respect to the treatment of potential backgrounds to microlensing. 
We note that an important source of background are supernovae in galaxies behind the 
LMC. 

Subject headings: dark matter — Galaxy: structure, halo — gravitational lensing — 
Stars: low-mass, brown dwarfs, white dwarfs 
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1. Introduction 

Following the suggestion of Paczynski (1986), sev- 
eral groups are now engaged in searches for dark mat- 
ter in the form of massive compact halo objects (MA- 
CHOs) using gravitational microlensing, and many 
candidate microlensing events have been reported. 
Reviews of microlensing in this context are given by 
Paczynski (1996) and Roulet & MoUerach (1996). 



Previously ( Alcock et al. 1997a ) we conducted 
an analysis of 2.1 years of photometry of 8.5 million 
stars, and found 6-8 microlensing events, implying 
an optical depth towards the LMC of 2.9j;n Q x 10~^ 



for the 8-event sample and 2.1 



+1.1 

-n 7 



-0.9 

10-^ for the 6- 



event sample (( Alcock et al. 1996a ; Alcock et al 



1997a ), hereafter A96 & A97, respectively). Inter- 
preted as evidence for a MACHO contribution to the 
Milky Way dark halo, this implied a MACHO mass 
out to 50 kpc of 2^1'^ X 10" M©. Depending on the 
halo model this meant a MACHO halo fraction of be- 
tween 15% and 1 00%, and a typical MACHO mass 
of 0.1 to 1 Mq ([Alcock et al. 1997a|; [Gates et al 



1995). The EROS group has reported 2 candidates 
(Aubourg et al. 1993; Ansari et al. 199^ , consis- 



tent with the above results. Recently the EROS's 
updated and expanded survey, EROSII, has reported 
two new events that they interpret as limiting the 
amount of halo dark matter (Milsztajn 1999), but are 
consistent with both A97 and the results of this pa- 



per. The OGLE ( [Udalski et al. 1997[ ) collaboration 
also reported one LMC microlensing event in 1999. 
All claimed LMC events have characteristic timescales 
between t ~ 34 and 230 days, while searches for short- 
timescale events with timescales 1 hour <^t<10 days 



have revealed no candida t es to date (Aubourg et al. 



1995[ ; [Alcock et al. 1996d[ ; [Alcock et al. 1998a[) allow- 



ing important limits to be set on low mass dark mat- 
ter. In addition, two candidates have been observed 



towards the SMC (Alcock et al. 1997b; 


Alcock et al. 


1998b; 


Albrow et al. 1999; 


Palanque-Delabrouille et 


al. 1998; 


Afonso et al. 1999|; ?; Udalski et al. 199^; 


Afonso et al. 1999 


), but the small number of events. 



location of the lenses, and large expected SMC self 
lensing rate reduce their usefulness as a probe of the 
dark halo. 

Conclusions based on our previous work suffered 
due to Poisson error as a result of the small number 
of events. Increasing the time span monitored from 
2.1 years to 5.7 years and increasing the number of 
monitored fields from 22 to 30, gives 13 to 17 events 



(depending on the cuts used), as well as greatly in- 
creases our sensitivity to long duration events and 
therefore to higher mass MACHOs. It also increases 
coverage over the face of the LMC, providing a use- 
ful tool to distinguish between various interpretations 
of the microlensing events. In addition, while our 
previous analyses contained the most careful evalua- 
tions of microlensing detection efficiency ever done, 
we have made several important improvements, and 
have thoroughly tested the robustness of our methods. 
In addition we have more fully and carefully investi- 
gated sources of potential background to microlens- 
ing, in particular LMC variable star background (i.e., 
bumpers A97) and supernovae in galaxies behind the 
LMC. For example, one candidate event classified in 
A97 as microlensing (LMC-10) is now removed as a 
probable background supernova. 

The nature of microlensing implies that many of 
our events will be low signal-to-noise, so to test the ro- 
bustness of our results and to estimate systematic er- 
ror due to our event selection methodology, we present 
two independently derived sets of selection criteria, 
with two corresponding sets of events and efficiency 
determinations. One set is designed to select only 
high signal-to-noise events and is modeled on the se- 
lection criteria used in A97. The other is designed to 
be inclusive of lower signal-to-noise events and also 
exotic microlensing events, and makes heavier use of 
several new statistics. While the number of microlens- 
ing candidate events selected by the two sets of cuts 
differ, the corresponding efficiencies compensate, and 
the resulting optical depth values, halo fractions, etc. 
are essentially the same. This suggests that the sys- 
tematic error in our optical depth, etc. due to our 
choice of cuts is small. Finally, we implement an 
improved likelihood analysis that self-consistently in- 
corporates currently available information on known 
stellar backgrounds. 

The increase in number of events, improved effi- 
ciency determination, and more thorough investiga- 
tion of systematic errors and backgrounds, such as 
bumpers and supernovae, makes the results of this 
paper the most accurate to date. At this point, uncer- 
tainties in the model of the Milky Way and the model 
of the LMC dominate both the quantitative and in- 
terpretational aspects of microlensing as a probe of 
dark matter. 

The plan of the paper is as follows: in § || we out- 
line the observations and photometric reductions. In 
§ 1^ we describe our microlensing event selection crite- 
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ria, present the resulting candidates, and discuss sev- 
eral sources of background to microlensing, including 
bumpers and supernovae. In § ^ we estimate our de- 
tection efficiency, which has been improved in a num- 
ber of ways. In § H we show the distributions of the 
selected events in the color-magnitude diagram, loca- 
tion on the sky, and impact parameter. We compare 
with predicted distributions, thereby testing the mi- 
crolensing hypothesis. In § ^ we provide various anal- 
yses of the sample. We calculate the optical depth, 
and discuss why it is a factor of two smaller than 
in A97. We perform a likelihood analysis that ex- 
plicitly includes models of the Milky Way and LMC 
stellar populations to find new estimates of the MA- 
CHO contribution to the dark halo and the number of 
expected events from known stellar populations. We 
find new favorable mass ranges for the lenses if they 
are halo objects. We also discuss various interpreta- 
tions of our results, including the possibility that no 
MACHOs exist and all the lensing is due to stellar 
lenses. 

Note that many of the reduction and analysis pro- 
cedures used here are very similar to those in A96 
& A97, to which we refer extensively. A more rig- 
orous description of our detection efficiency, which is 
only briefly outlined in this paper, may be found in 
the companion paper, Alcock et al. (2000a), and in 
Vandehei (2000). The reader is encouraged to consult 
these papers to understand the details of the exper- 
iment, but we will repeat the main points here for 
clarity. 

2. Observations and Photometric Reductions 

The MACHO Project has had full-time use of the 
1.27-meter telescope at Mount Stromlo Observatory, 
Australia, since 1992 July. Observations are sched- 
uled to be completed at the end of 1999 December. 
Details of the telescope system are given by Hart et 
al. (1996) and of the camera system by Stubbs et al. 
(1993) and MarshaU et al. (1994). Briefly, corrective 
optics and a dichroic are used to give simultaneous 
imaging of a 42 x 42 arcmin^ field in two colors, us- 
ing eight 2048^ pixel CCDs. As of 1998 March, over 
70,000 exposures had been taken with the system, 
over 5 TBytes of raw image data. About 55% are of 
the LMC, the rest are of fields in the Galactic center 
and SMC. 

In this paper, we consider the first 5.7 years of data 
from 30 well-sampled fields, located in the central 5° x 



3° of the LMC; field centers are listed in Table ^ and 
shown in Figure |^. 

The observations described here comprise 21,570 
images distributed over the 30 fields. These include 
most of our observations of these fields in the time 
span of 2067 days from 1992 September 18 to 1998 
March 17 as well as a fraction of our observations 
taken between 1992 July 22 and 1992 August 23 when 
our system was still in an engineering phase. The 
mean number of exposures per field is 21570/30 = 
719, with a range from 180 to 1338. The sampling 
varies between fields (Table p, since the higher pri- 
ority fields were often observed twice per night with 
an average of about 4 hours between exposures. 

The photometric reduction procedure was very 
similar to that described in A96 & A97; briefly, a 
good-quality image of each field is chosen as a tem- 
plate and used to generate a hst of stellar positions 
and magnitudes. The templates are used to "warm- 
start" all subsequent photometric reductions, and for 
each star we record information on the flux, an error 
estimate, the object type, the of the Point Spread 
Function (PSF) flt, a crowding parameter, a local sky 
level, and the fraction of the star's flux rejected due 
to bad pixels and cosmic rays. Details of the MACHO 
image data and photometry code (SoDOPHOT) are 
provided in Alcock et al. (1999c). The resulting data 
are reorganized into lightcurves, and searched for vari- 
able stars and microlensing events. The LMC 5.7-year 
photometry database is about 200 Gbytes in size. 

For the 22 flelds reported on in A97, we have well 



calibrated photometry (Alcock et al. 1999c), but for 
the eight new flelds our photometry has been only 
roughly calibrated on a global basis (see § 2 in A97). 
Event selection is generally based on this rough cali- 
bration, but as noted below, we report the well cali- 
brated magnitudes and colors when possible. 

We have corrected a minor complication in A97 
where, for software-related reasons, we had used dif- 
ferent templates for the flrst and second year's reduc- 
tions of 6 of our fleldsf^ For these fields there was not 
a one-to-one correspondence between the set of stars 
in the 2 distinct years used, and the first and second 
years had to be analyzed separately. All photometry 
in these 6 fields have been re-run using the new gener- 
ation of templates and the lightcurves have now been 
merged onto a common photometric system. 



''Fields 1, 7, 9, 77, 78, & 79 
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3. Event Detection 

The data set used here consists of about 256 biUion 
individual photometric measurements. Discriminat- 
ing genuine microlensing from stellar variability, back- 
ground, and systematic photometry errors is hard, 
and the significance of the results depends upon the 
event selection criteria. 

The selection criteria should accept 'true' microlens- 
ing events, and reject events due to intrinsic stellar 
variability and instrumental effects. The determina- 
tion of our event selection criteria could not be made 
before looking in detail at the lightcurves. We had to 
discover various background sources and learn how to 
perform event selection from the data we gather itself, 
making the selection criteria dependent on the data. 
As much as possible we have tried to base the selec- 
tion criteria upon our Monte Carlo artificial events 
(see § ^) . Although this allows us to place cuts along 
natural breaks in parameter space (which lessens the 
sensitivity of the final results on the exact placement 
of the cut) it did not allow us to fully explore the 
background of variable stars. This adds some sub- 
jectivity to our analysis, which we quantify below by 
considering two limiting cases. 

For each lightcurve, we compute a set of over 
150 temporal variability statistics. We use two lev- 
els of statistics: level-1 statistics are calculated for 
all stars, while level-2 statistics are calculated only 
for those stars that pass the level-1 selection criteria. 
We have developed selection criteria ( "cuts" ) that use 
the level-2 statistics to distinguish microlensing from 
backgrounds such as variable stars and noise. The 
selection criteria have evolved over the course of the 
experiment. As the volume of data on a lightcurve 
increases, the meaning of some statistics change in 
subtle ways. Thus one must be careful not to blindly 
apply selection criteria from one data set to another. 
For example, a fit ^ to a constant flux star over 2.1 
years of data will not necessarily be the same when 
computed using 5.7 years of data, due to changes in 
weather patterns (and thus seeing & sky level) and 
the CCD camera over the course of the experiment. 
In addition, in order to increase our sensitivity to low 
signal-to-noise and exotic microlensing, our level 1 cri- 
teria (see below) have been loosened relative to those 
used in A97. This means that our set of level-1 can- 
didates contains more variable stars and other noisy 
events, thus requiring changes to the final level-2 se- 
lection criteria. Because of the changes in the level- 



1 criteria, the selection criteria used in A97 are no 
longer appropriate for the year 5.7 data. About 45 
lightcurves, 26 of which are clearly noise or variable 
stars, would pass the A97 criteria applied to the cur- 
rent data. 

We select the events using two different and inde- 
pendently developed sets of level-2 selection criteria. 
This allows us to explore possible systematic error 
due to the choice of cuts. While the goal of both sets 
of cuts is to select as many microlensing lightcurves 
as possible, while rejecting as many non-microlensing 
lightcurves as possible, the two sets of cuts were ex- 
plicitly developed with complementary philosophies 
in mind. The first selection criteria (hereafter referred 
to as criteria A) was designed to be rather tight, only 
accepting events with a single highly significant bump 
in either passband, while requiring the baseline to re- 
main very flat, as expected in simple microlensing. 
These cuts resemble those in A97 and for the most 
part statistics similar to those described in A97 were 
employed. The second selection criteria (criteria B) 
was designed to be rather loose, in an attempt to 
search for exotic or low signal-to-noise microlensing 
candidates. This second set of cuts also looked for 
a single significant bump in cither passband with a 
fiat baseline, but made use of some new statistics not 
available in A97. The new statistics (described below) 
better characterize and filter out some variable stars 
and noisy events. To conservatively estimate the sub- 
jective nature of our event selection, marginal events 
suspected of being supernova, etc. are preferentially 
rejected from set A, but kept in set B. 

Note that as long as the experiment's event detec- 
tion efficiency is calculated properly, and the selec- 
tion criteria are sufficiently stringent to accept only 
real microlensing events, changes in the selection cri- 
teria should be accounted for in the efficiency calcu- 
lations, and the details should not greatly affect the 
final results (in the limit of a large number of detected 
events). This statement implicitly ignores exotic mi- 
crolensing events such as binary lens or parallax mi- 
crolensing events which have lightcurves that differ 
from those used in our efficiency determination. We 
have not determined our efficiency for exotic events, 
but selection criteria B is designed, in part, to be more 
sensitive to such events. Furthermore, a much more 
sensitive search for exotic lensing events has been car- 
ried out, and none were found. So, we do expect that 
the difference between selection criteria A & B is a 
reasonable indication of our selection criteria system- 
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atic error, and we find that this difference is fairly 
small as outlined here and discussed in § ^. 

We have summarized the old (A97) and new selec- 
tion criteria (A & B) in Table |[ In § 3.1 we briefly 
describe some of the analysis and statistics used by 
these selection criteria and in § we present the set 
of events selected by criteria A & B. In § |3.3| we iden- 
tify two main sources of background to microlensing 
and discuss how they can be removed from the true 
microlensing pool. Finally, in § 3.4 we remove our 
identified background events and summarize the final 
sets of microlensing candidates (the final set A & B) 
used to compute the results of this paper. 

3.1. Selection Criteria Statistics 

Photometric measurements with questionable PSF 
fit, too much crowding, missing pixels, or cosmic rays 
are flagged as suspect and removed from further con- 
sideration. The event detection then proceeds in two 
stages. The first stage, defining a level-1 collection of 
candidate events, is similar to that described in A97; 
a set of matched filters of timescales 7, 15, 45, and 100 
days is run over each lightcurve. If after convolution, 
a lightcurve shows a peak above a pre-defined signif- 
icance level in either color, it is defined as a 'level-1 
candidate'. We also make use of a new filter that looks 
for bumps of any duration and add these lightcurves 
to the level-1 pool of candidates. (We found no addi- 
tional candidates by this change, however.) For level- 
1 candidates, a full 5-parameter fit to microlensing is 
made, and many level-2 statistics describing the sig- 
nificance of the deviation, goodness of fit, etc. are cal- 
culated. We use the standard point-source, point-lens 



approximation (Refsdal 1964; Alcock et al. 1996a). 
The 5 free parameters of the fit are the baseline flux 
in red and blue passbands foR, foB, and the 3 param- 
eters of the microlensing event: the minimum impact 



parameter in units of the Einstein radius, 



the 



Einstein diameter crossing time, t = 2rE/vj_, and the 
time of maximum magnification, ^max- Later, instead 
of Mmin wc will ofteu usc the fit maximum magnifica- 
tion A,„ax = ^(wmin), which is more closely related 
to the observed light curve. 

Lightcurves passing loose cuts on these statistics 
are defined as 'level-1. 5' candidates, and are output 
as individual files along with their associated statis- 
tics. In the present analysis, there are approximately 
150,000 level-1. 5 candidates, almost all of which are 
variable stars or noise. 



Some of the important statistics used by both 
criteria A & B are the chi-squares of various fits. 
For example, a powerful signal-to-noise statistic is 
Ax^ = x?onst - xLi where Xconst and xLi are the 

values for the constant flux and microlensing fits, 
respectively. Ax^ is the effective 'significance' of the 
event summed over all data points. Xpg^k refers to 
the of the microlensing fit in the "peak" region 
where Afa > 1.1. A reduced x^ of microlensing 
fit outside the interval i,nax i 2t is also computed, 
^mi-out/'^dof- Other useful statistics include the av- 
erage values of crowding for a star, the microlensing 
fit values, the magnitude and color, the number of 2- 
sigma high points in the peak region, the number of 
points on the rising and falling side of the peak region, 
and the number of points outside the peak region. 

Based upon our experience gained in A96 & A97 
we have developed a number of new statistics. One 
such statistic is the fraction of points in the peak that 
lie above the lightcurve's median, Mii/iVpk. This new 
statistic is useful for removing events with spurious 
deviant points associated with crowding/seeing in- 
duced fluctuations, satellite/asteroid tracks and other 
causes. To further help in removing similar spuri- 
ous events we also compute the fraction of points 
in the peak rejected due to bad PSF measurements, 
pkpsfrej, and large crowding values, pkcrdrej. 

We have found it beneflcial to concentrate on 
statistics that help in rejecting variable star back- 
ground. One such statistic is the ratio of powers 
in the two passbands, bauto/rauto. The quantity 
bauto is the sum of the absolute values of the flux 
away from the median in the blue bandpass, and like- 
wise for rauto for the red bandpass. Our variable 
star background generally has more power in the red 
than in the blue, while our Monte Carlo microlensing 
events tend to have equal power in both passbands, 
even when heavily blended (that this is the case is 
not surprising given the fact that the color difference 
between any two stars in the LMC is somewhat re- 
stricted) . A second new statistic on the ratio of pow- 
ers in each filter, pfwsr is also computed, but uses 
only points measured simultaneously in both filters 
and is thus normalized differently than bauto/rauto. 
To further help reject variable stars we compute a 
cross-correlation coefficient between the red & blue 
filters, rbcrossout, but use only data points that are 
outside the peak of the event, ±1.0i, to avoid the 
peak contributing. A powerful new statistic for re- 
jecting variable stars is a robust reduced chi-square 
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fit to a constant flux, Xrobust-out/^dot j also computed 
outside the interval ±1.0t . 

Another new signal-to-noise statistic we compute 
is pfrdev, which is a from the baseline counting 
only upward excursions in the filter window and sub- 
tracting a penalty per measurement in the bin so that 
2-sigma points break even. We subtract the largest 
single date contribution (1 or 2 points) to ensure ro- 
bustness against single strongly deviant observations. 
We also compute a similar for the second most sig- 
nificant non-overlapping filter trigger, pfrdev2. The 
statistic pfrdev is very similar to Ax^/(Xmi/-^dof ); 
and a comparison of the two on both real data and ar- 
tificial data yields similar results, even though the lat- 
ter is derived using the microlensing shape. This gives 
us added confidence that a cut on Ax^/(Xmi/-^dof) 
is not very 'shape' dependent. We have opted to use 
^X^/ (Xmi/-^dof) throughout, but make use of the sec- 
ond peak significant pfrdev2 to ensure the unique- 
ness of the event. This statistic is useful for elimi- 
nating variable stars, but must be used with caution 
to avoid missing exotic lensing events such as binary 
lenses or sources that could exhibit a second 'bump'. 
Similarly, cuts on Xmi-out/^dof and xLbust-out/^dof 
could potentially bias against detecting widely sepa- 
rated binary microlenses. We ran our selection cri- 
teria with and without these cuts to ensure that no 
exotic lensing events were missed. No additional can- 
didates were found. 

Our set A selection criteria are designed to ac- 
cept high quality microlensing candidates, while us- 
ing mainly the statistics described in A97. Some of 
the cuts on these statistics have been loosened as we 
have developed a better understanding of our variable 
star background; for example the cut on magnification 
was loosened from ^max > 1-75 to its present value 
^max > 1.49, because our main background of vari- 
able stars, a class of blue variables called bumpers (see 
below), almost never show fit magnifications larger 
than 1.5, and are well isolated in the color magnitude 
diagram. We also loosened our main significance cut 
AxV(xii/^dof) < 400 (from 500 in A97) because a 
number of other statistics were tightened. These in- 
clude, Ainax > 3^ (from 2a where a is the average 
red/blue error in magnitudes) which is a signal-to- 
noise cut, Xmi-out/^dof < 1-8 to reject variables (from 
< 4.0 in A97) and AxV(Xpoak/^dof) > 350 (from 
> 200 in A97), another signal-to-noise cut. This later 
cut is somewhat more reliant on the shape of the 5 pa- 
rameter microlensing fit and increases the likelihood 



of rejecting exotic microlensing; for this reason, crite- 
ria B below does not use this cut. 

Our set B selection criteria are designed to accept 
any lightcurves with a significant unique peak and a 
fairly fiat baseline. Our selection criteria B are sum- 
marized in Figure ^ which illustrates our two most 
important signal-to-noise cuts: the cut on magnifica- 
tion Amax and the cut on AxV(Xpoak/^dof)- Events 
that passed the basic cuts (all cuts used by criteria 
B minus the ordinate & abscissa cuts of Figure ^ 
are shown as solid dots and labeled. The final cut on 
magnification A,„ax and Ax^/(Xpcak/^dof ) are shown 
as solid lines (for comparison dotted lines for criteria 
A are also illustrated). Open circles indicate events 
which fail criteria B's bumper cut and solid squares 
events which fail criteria B's uniqueness cut. The 
new statistics on the number of points rejected in the 
peak region, pkpsfrej +pkcrdrej, and the fraction 
of points above the baseline in the peak, A^hi/A^pk are 
useful for eliminating spurious noise- induced events. 
New statistics on uniqueness, pfrdev2, the passband 
power ratio pfwsr, the red/blue cross-correlation co- 
efficient outside the peak rbcrossout, and the ro- 
bust Xrobust-outZ-^dof baselinc statistic are useful in 
removing periodic and quasi-periodic variable stars. 
We have also more carefully characterized our main 
source of variable star background, the bumpers, in a 
magnitude-color-magnification space (see Figure 
With the background of variable stars more effec- 
tively removed we can reduce both our reliance on 
any 'shape' dependent criteria and lower the signifi- 
cance level of a detection. In § || we demonstrate the 
relative looseness of criteria B over A as well as the de- 
creased dependence on shape. However, one potential 
difficulty with criteria B is its inability to discriminate 
against some types of variable stars, such as CVs and 
SNs, that might exhibit strongly asymmetric and/or 
chromatic lightcurves during the 'event' but remain 
constant for long periods of time. Supernova removal 
is discussed separately in more detail below. See Ta- 
ble |2| for a complete list of the individual cuts used by 
selection criteria A and B. 

3.2. Microlensing Candidates 

We find 19 lightcurves that pass criteria A and 
29 lightcurves that pass criteria B (before apply- 
ing the supernova cuts described below). All the 
lightcurves passing criteria A also pass criteria B. The 
29 lightcurves are shown in Figure ^ and their mi- 
crolensing fit parameters are listed in Table ^. Events 
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that do not pass criteria A are marked with an aster- 
isk. Parameters for fits including the possibihty of 
blending with an unlensed star in the same seeing 
disk as the lensed star are given in Table g The un- 
blended fits are displayed as a thick line in Figure ^, 
used for all statistics except comparison with super- 
nova. Note that our events here are numbered as in 
A97 to avoid any possible ambiguity. Thus the first 
event described here that was not described in A97 
is event 13. Finding charts for the events, as well as 
the full lightcurves can be found on the internet at 
http : / / wwwmacho . anu . edu . au/ 

Six of these 29 lightcurves (la, lb, 10a, 10b, 12a 
and 12b) actually correspond to only three stars that 
occur in field overlap regions; the two lightcurves for 
each star are based on independent data and reduc- 
tions. Two lightcurves are also of the same star (7a 
and 7b) but were not in field overlaps. Event 7 was 
bright enough and in a locally crowded enough region 
that some of the flux from the primary (7a) contam- 
inated a secondary (7b) neighbor causing a spurious 
detection. Event 7a passed both criteria A & B, while 
event 7b passed only criteria B due to its low signal- 
to-noise. Thus before supernova removal there are 16 
unique events found by the criteria A, and 25 unique 
events found by the criteria B. 

3.3. Background 

3.3.1. Bumpers 

As noted in A97 a potential source of background 
to microlensing is a class of bright blue variables 
which we refer to as bumpers. Although associated 
with Be stars, which are known to show periodic out- 
bursts in our Galaxy, the true nature of these variable 
stars is still unknown. However, it is possible to elim- 
inate bumpers as a serious source of background as 
they can be well isolated in a multi-parameter space. 
For example, microlensing fits to bumpers seen in our 
data almost never return magnifications larger than 
1.5, typically much less, as can be seen in Figure ^ 
(open circles). 

One does not have to be confined to magnifica- 
tions above 1.5 as the bumpers are also well isolated 
in the color-magnitude diagram (CMD) as illustrated 
in Figure ^. Here a typical CMD of the LMC is shown 
with a scattering of small dots. The 29 lightcurves 



'A mirror site exists at 

http://wwwmacho.iiicmaster.ca/ and the site for microlensing 
alerts is http://darkstar.astro.washington.edu. 



that pass selection criteria B are shown as filled cir- 
cles and labeled. As in Figure ^ open circles indi- 
cate events which fail criteria B's bumper cut. Filled 
boxes indicate events which fail criteria B's magni- 
fication cut ^max > 1-34. The final cut on bright- 
ness V > 17 and color-magnitude-magnification (the 
'bumper' cut) are marked as solid lines for crite- 
ria B (for comparison the dotted line illustrates the 
'bumper' cut for criteria A). If a potential event falls 
within the boxed region labeled 'Anax > 1-75' then 
it must have a magnification greater then 1.75 to be 
included in the criteria B set of events. We have visu- 
ally inspected the ~ 300 bumper candidate lightcurves 
and confirm that the vast majority of them show the 
slight asymmetry and other characteristics typical of 
the bumpers described in A96. In fact, due to the cut 
on uniqueness these bumpers have only single bumps 
and evidently represent one-time only bumpers or 
bumpers with inter-bump intervals longer than 5.7 
years. 

3.3.2. Supernovae 

Another serious source of potential contamination 
in microlensing surveys, that has not been given suffi- 
cient attention before, are supernovae (SN) occurring 
in galaxies behind the LMC. These background SN 
are picked up in the crowded fields, and their host 
galaxies are not always easy to identify in ground- 
based images. The fact that they occur only once 
and show a flat baseline before and after the 'event' 
make SN interlopers a serious concern. 

A first step in understanding this source of con- 
tamination is to estimate the number of SN we might 
see during the course of the experiment. Because of 
the recent interest in SN type la as standard candles, 
the rate of SN (both type I and II) occurring in field 
galaxies is now fairly well known. We use a typical 
rate of 0.5 SN /year/D° with peak m agnitude brighter 
than V ^20 ( [Woods fc Loeb 1998]) . The duration of 
the experiment is 5.7 years and covers 13.5 which 
suggests we should have approximately 38 SN in our 
data set. This does not include our SN detection ef- 
ficiency and so is an overestimate. We expect our 
efficiency for detecting SN to be on the order of 5- 
15%, similar to that of detecting microlensing events 
(see § ^) , due to the similar shapes of the correspond- 
ing lightcurves, implying we are likely to see ^2 — 6 
SN in the current data set. 

If the density of galaxies behind the LMC is aver- 
age, then the probability of finding a galaxy in prox- 
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imity of a star should be low. Therefore, a robust 
way of eliminating potential SN interlopers would be 
a search for a background host galaxy in an image. 
If such a galaxy were found within a pre-determined 
radius, the microlensing interpretation would be un- 
likely. Since the area covered by nearby galaxies is 
small, the correction to our experimental efficiency 
would be small. Unfortunately, our ground-based im- 
ages do not have the resolution and low enough sky 
count to perform such a search with confidence. We 
do have HST observations of 8 of our events (see Ta- 
ble ^ and these give us high confidence that 6 of them 
are not SN as there are no obvious background galax- 
ies anywhere near these events. Two events for which 
we have an HST image (events 10 and 12) do show a 
fairly bright spiral galaxy within 5 arcsecs. In retro- 
spect, the host galaxy for event 10 shows up in our 
ground-based images but could not have been recog- 
nized as a galaxy without further data. This event 
was noticeably asymmetric in A97 and was classified 
as a marginal microlensing candidate. In calculating 
the optical depth, it was rejected from the 6 event 
sample, but included in the 8 event sample. Although 
the effect of its inclusion was small, increasing the re- 
ported optical depth of the 8 event sample by only 
7%, its presence underscores the need to take the po- 
tential SN contamination seriously. 

Unfortunately, with an incomplete sample of HST 
images we are unable to follow the above prescrip- 
tion to reject SN interlopers completely. We have 
requested HST time to observe the remaining events, 
and have examined the best available ground based 
images for evidence of a background galaxy. How- 
ever, we can also use the shape of the lightcurve to 
help distinguish SN from microlensing. Even though 
the lightcurvcs of SN type II are not well understood 
and exhibit a wide range of behavior, SN of Type la 
have been studied in detail. Their lightcurves are very 
similar once distance , reddening, and a shape param - 
eter are allowed for ( Phillips 1992 ; Riess et al. 1995 ) . 
In addition, type la's pr obably dominate the SN rate 
in flux limited samples ( Woods fc Loeb 1998 ) . Using 
the SN type la templates of Riess, Press, & Kirshner 
(1996) transformed to the MACHO photometric sys- 
tem, we have applied a 6 parameter SN type la fit to 
all 29 lightcurves. The 6 free parameters of the fit are 
the baseline flux of the photometered object in red 
and blue passbands /qa, fovi and 4 parameters that 
describe the SN type la lightcurve: the time of peak 
ipcak, the distance modulus in red band /i^ and blue 



band /.ty (fit independently to account for the possi- 
bility of reddening) and a shape parameter 5 which 
parameterizes how SN of type la become longer in 
duration when intrinsically brighter. Emperically, we 
discovered that by allowing any acceptable value of 5, 
in many cases our best fit SN shape was well outside 
the range of observed SN (e.g. (5 = 5). This was espe- 
cially true for high-quality microlensing events, where 
the SN template provides a poor fit to the shape of 
the lightcurve. We therefore limit our range in 5 to 
be between -0.5 and 0.75, and note that for most of 
the events with poor SN fits, we find 5 pegged at one 
of these values. The SN type la fits are displayed as 
a dashed line for the events we categorize (see below) 
as SN in Figure H] and the fit parameters for all 29 
lightcurves are given in Table |. 

We summarize all the relevant available informa- 
tion for the 29 events in Table |[ Columns 2-4 com- 
pare the blended microlensing fits XMb/^dof with 
the SN type la fits XsNZ-^dof- A positive value of 
^XsN-ML indicates a better fit to blended microlens- 
ing, while a negative value a better fit to SN type la. 
Inspection of Table || reveals 10 lightcurves (8 events) 
that are better fit by SN type la. As a consistency 
check column 5 indicates the presence or absence of an 
obvious background galaxy within 10 — 15 arcsecs 
as determined using the best available image of the 
event (image source given). This background galaxy 
identification is subjective. Part of the problem is 
the severe crowding of the ground-based images used 
(MACHO or CTIO images). In every case, where 
a galaxy is probably present the fit to a SN type 
la template is preferred over the blended fit, giving 
us some assurances of the overall correctness of the 
procedure. We have no spectra to confirm the hy- 
potheses that any of these 8 events are, in fact, SN 
of type la. But given the fact that we should see 2-6 
SN in our survey we feel it is nevertheless prudent to 
eliminate these 8 events as potential interlopers. We 
thus implement as our final cut (on both selection 
criteria A & B) the requirement that the blended mi- 
crolensing fit be preferred over the SN type la fit, or 
AxsN_ML > O-O- This eliminates events 10 and 12 
from set A and events 10, 11, 12, 16, 17, 19, 24, and 
26 from set B. Event 22, clearly not a type la, is a 
special case and will be discussed later. The advan- 
tage of using this simple cut is that we may quantify 
the effect it has on the detection efhciency (see § 
The effect is negligible since less then 0.7% of artifi- 
cial standard microlensing events are falsely rejected 
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by this cut. 

There are several important caveats to using the 
SN la fits to reject potential SN interlopers. First. 



hand, it is also possible that event 22 belongs to a 
class of SN ("slow" type I l n) similar to SN 1988Z 
(iRosecherr k Schaefer 1999| ; ^chlcgel 1990| ; ^tathakis 



both microlcnsing and SN come in a varietv of fla- 



vors. Although type la have well defined lightcurves, 
other types of SN are not so well behaved (types lb, 
Ic, Hp, III, etc.) and can come in a variety of dura- 
tions and asymmetries. We note here that two of our 
SN interlopers (events 11 and 24) are likely to be SN 
type Hp as judged by the plateau seen in both pass- 
bands about 25-50 days after maximum (never seen 
in both passbands for type la). However, even these 
lightcurves are usually better fit by SN type la than 
blended microlensing, mostly due to the asymmetry 
that SN lightcurves typically exhibit. 

Second, exotic microlensing such as found in binary 
lenses or parallax events, could mimic the asymmetry 
of type la SN. There are several reasons why we be- 
lieve that this is not a major problem: 1) event 9, 
a binary lens event, is better fit by a single lens mi- 
crolensing lightcurve than a SN la lightcurve, 2) out of 
many microlensing events towards the Gal actic bulge 
les s than 10% are of clearly exotic type ( Alcock ei 



al. J_999b)J^, 3) exotic microlensing should show the 
"wrong" sign of asymmetry 50% of the time, and we 
have no examples of this among our events. 

However, two events do stand out as potentially 
worrisome. Event 26 is better fit by SN type la, but 
does not show a clear background galaxy in our deep- 
est CTIO 0.9m images. It could be an example of 
exotic lensing. For purposes of this paper, we reject 
this as a potential SN interloper, but must await a 
better image before making any definitive conclusions 
on this event. Event 22 is our longest duration event 
and is clearly asymmetric to the eye. However, even 
with the clear asymmetry this event is better fit with 
blended microlensing than with a SN type la, mainly 
because type la's are not observed to last this long. 
The asymmetry of this event is well fit by microlensing 
parallax, which would be a natural explanation given 
the very long duration of this event. There is also no 
strong evidence of a galaxy in ground-based images, 
though the object is slightly extended. Microlensing 
parallax is a good explanation for this event, and if 
this interpretation is correct, it appears that the lens 
is probably a w hite dwarf in a flattene d halo or thick 
disk population( Alcock et al. 2000b| ) On the other 



Sadler 1991). Therefore, in the spirit of our selec- 



^^However, the lens populations towards the LMC may be differ- 
ent than that towards the bulge. 



tion criteria, we reject event 22 from the exclusive set 
A and keep it in the inclusive set B. Since event 22 is 
our longest duration event and therefore contributes 
maximally to the optical depth, this is also the con- 
servative approach to exploring the sensitivity of our 
results to the selection criteria. 

3.4. Hovif Many Events? 

Here we briefly summarize the events that will con- 
stitute set A & set B. The main results of this paper 
rest upon these two sets of candidate microlensing 
events. 

Criteria A selected 19 lightcurves corresponding to 
16 unique events (events 1, 10, & 12 in field overlaps). 
Of these 16 unique events two (events 10 & 12) are 
rejected as SN interlopers due to the final SN cut 
AxsN-ML ^ ^^'^ unexplained event 22 is re- 
jected in the spirit of criteria A being exclusive. This 
leaves set A containing 13 events: 1, 4-8, 13-15, 18, 
21, 23, & 25. 

Criteria B selected 29 lightcurves corresponding to 
25 unique events (events 1, 10, & 12 in field overlaps 
and event 7 duplicated via contamination). Of these 
25 unique events 8 (events 10, 11, 12, 16, 17, 19, 24, 
& 26) are rejected as SN interlopers. Here, we do not 
reject event 22 but leave it in the set as a potential 
exotic lensing event. As a result, set B contains 17 
events: 1, 4-9, 13-15, 18, 20-23, 25, & 27. 

4. Detection Efficiency 

The detection probability for individual events de- 
pends on many factors, e.g. the 3 event parameters 
^max, t , imax, and the unlensed stellar magnitude, 
as well as our observing strategy and weather condi- 
tions. Such a complicated dependence is most natu- 
rally found with a Monte Carlo technique. We may 
simplify the dependence to some extent by averag- 
ing over the known distributions in Amax, ^max, the 
stellar magnitudes, and the known time-sampling and 
weather conditions, to derive our efficiency as a func- 
tion only of event timescale, £{t). 

We have computed our detection efiiciency using 
a method similar to that outlined in A97, but with 
a number of improvements. A full discussion of the 
method, with detailed results, is given in Alcock et 
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al. (2000a) and Vandehei (2000). Briefly, we generate 
simulated microlensing events with t logarithmically 
distributed in the range 1-2000 days over the slightly 
wider time interval, (JD - 2,448,623.5) = 190.0 to 
2277.0, and add these simulated events into the ex- 



ten led time span of observations. A large database 



of artificial star tests is used to simulate the effects 
of blending. A number of systematic photometric ef- 
fects, including the response of flux, error bars, and 
the photometric flags outlined in § || are also included. 
The Monte-Carlo procedure takes into account the 
actual spacing and error bars of the observations, so 
any variations in sampling frequency, weather, seeing, 
etc. between the A97 data set and the current data 
set are automatically accounted for. 

One of the primary shortcomings of the efficiency 
analysis presented in A970 was a lack of faint 'stars' 
in the artificial star tests, which are used to add sim- 
ulated events. In the current analysis, a large num- 
ber of faint stars down to V ^ 24.5 are used. This 
is 2.5 magnitudes fainter than in A97. It is also 2.5 
magnitudes fainter than our faintest detected objects. 
The present analysis also makes use of a much larger 
database of artificial stars (x5) sampled over a larger 
(xl5) and more fairly distributed set of observing 
conditions (stellar density, seeing, & sky). Another 
major improvement, not fully recognized as a major 
source of uncertainty in A97, is a normalization of 
our fields to the underlying luminosity function. A 
main issue in efficiency determination is the distinc- 



tion between "objects" and stars. Objects are flux 



concentrations identified by the photometry code as 
stellar-like objects. Each object is typically a blend 
of many underlying LMC stars, any one of which can 
undergo microlensing, and it is important to identify 



the correct density of underlying LMC stars in each 



of our fields. 

The present efficiency is based on all stars in our 
fields, even those not uniquely identified because of 
signal-to-noise or crowding effects. These are ac- 
counted for by integrating the detection efhciency per 
star over the true underlying luminosity function (LF) 
in the LMC. Our LF in the LMC is derived from 
a combination of ground-based MACHO photometry 
(for stars with V < 20) and multiple HST WFPC2 
fields in the LMC bar (for stars with V > 20) and is 



^^We consider here only the 'photometric' efficiency defined in 
A97. The 'sampling' efficiency, also described in A97 is of little 
value for the discussion of this paper. 'Sampling' efficiencies for 
the present analysis may be found in Alcock et al. {2000a). 



described in Alcock et al. (2000a) and Alcock et al. 
(1999c). The shape of the LF appears universal in 
most of our fields and is well constrained for V < 22. 
But, any reasonable deviation from the adopted shape 
down to V ^ 24 has little effect on our efficiency (Al 



cock et al. 2000a) 



An important and yet uncertain factor is the nor- 
malization of the LF in each of our fields, which de- 
termines our effective sensitivity or exposure in star- 
years. Because our exposure in object-years is well 
known and the efficiency should converge at some 
magnitude (as fainter stars contribute less and less) 
we have chosen to factor this normalization into the 
efficiency. Thus our efficiency is properly defined 
as £{t) = {S/0){V < 24)£:stars(l^ < 24,f), where 
fstars(^ < 24, i) is the efficiency per star to V = 24 
and {S/0){V < 24) is the normalization, or the num- 
ber of stars per object (defined as the ratio of the 
true number of stars with < 24 to the number of 
SoDoPhot objects). The value of {S/0){V < 24) is 
10.84 ±2.4 stars per object and represents a weighted 
average over all 30 fields. The limiting magnitude 
V = 24: was chosen because the efficiency for du- 
rations t < 300 days converges. The efficiency for 
longer durations does not converge hy V — 24, and 
thus these are likely underestimated, with the under- 
estimation becoming worse for longer duration events. 
The exact point and speed of convergence is sensitive 
to the assumed shape of the LF and the cuts used, 
with criteria A converging more rapidly (see (Alcock 
et al. 2000a)). Since none of our 13-17 events have 



durations longer than 300 days our efficiency deter- 
mination for them is sound. The uncertainty in the 
efficiency is dominated by the uncertainty in the nor- 
malization, which we estimate to be ~ 20% (see (Al 



cock et al. 2000a) for a more detailed error budget). 



Efhciency results are shown in Figure |^. Selection 
criteria A is shown as a .solid line and criteria B as a 
dotted line. Also shown for comparison are the effi- 
ciencies presented in A96 {long dash) and A97 {short 
dash). Note that the efficiencies presented here have 
been scaled by a normalization term that accounts for 
the increased exposure due to all stars in our fields, 
down to y = 24 (as described above) . Strictly speak- 
ing, the efficiency defined in this manner is not con- 
strained to lie below one, though in practice it always 
does. This efficiency is defined relative to an 'expo- 
sure' of = 6.12 X 10^ object-years, which arises as 
follows: there are 11.9 million lightcurves in our total 
sample, and 20% occur in field overlaps. The relevant 
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time span is the 2087-day interval over which we add 
the simulated events; thus the exposure is 10.7 x 10^ 
objects x2087 days = 6.12 x 10'' object-years. This 
exposure is 3.4 times larger than in A97. Note that 
number of stars in field overlaps has increased from 
12% in A97, due to additional fields and a more care- 
ful calculation of the field overlap size. 

The most striking difference between the previous 
two data sets (A96 & A97) and the current set is the 
much higher efficiency at long durations. Most of this 
difference is a refiection of having 5.7 years of data 
instead of 2.1 years. Explicit cuts in A97 were made 
that removed any events with t > 300 days, while 
the current cuts both use t > 600 days. We also 
made a slight modification to our observing strategy 
intended to increase our sensitivity to long duration 
events. However, some of the difference also lies in a 
quirk of the year 2 data set. In A97, six of the dens- 
est fields had their lightcurves cut in half, roughly, 
due to an early generation of templates used to re- 
duce the photometry for these fields. As described 
in § H the photometry for these six fields has been 
re-run with the new generation of templates and the 
lightcurves merged onto a common photometric sys- 
tem. The effect of this was a lowered efficiency in A97 
for events with durations longer than t ^ 100 days 
(mainly due to the required 40 baseline points in the 
'halved' lightcurves). The problem did not exist in 
the year 1 data, thus the rather similar behavior of 
year 1 and year 2 at large t , even though the later 
had twice the coverage. 

The relative looseness of selection criteria B over 
criteria A discussed in § |3.l| (i.e. 17 vs. 13 events) 
is well illustrated in Figure 1^. Only for short dura- 
tions, t < 10 days, is criteria A more efficient, due 
to criteria B's larger number of required high points 
(> 10) as compared with criteria A (> 7). Less than 
half of the difference in efficiency between criteria A 
& B is explained by the different Amax cuts; given 
Anax cuts of 1.49 & 1.34 for criteria A & B respec- 
tively, we naively expect criteria B to recover 17% 
more events. In fact, only event 27 (7%) would have 
been missed had criteria B's vlmax cut been increased 
to 1.49. The remainder of the difference lies primar- 
ily in the effect of two cuts. Many of the events that 
failed criteria A did so because they failed either the 
cut in AxV(Xpoak/^dof) or the cut in Xmi-out/^dof • 
Both of these cuts have been tightened from their 
year 2 values. From our Monte Carlo events we 
find that tightening these two cuts has the follow- 



ing two effects: (1) lessens our sensitivity to moder- 
ately or strongly blended events, and (2) tightening 
the Ax^/(Xpoak/^ciof ) cut also decreases sensitivity to 
exotic microlensing and other asymmetric lightcurves, 
such as supernovae. For example, the slightly asym- 
metric event 26 did not pass criteria A because of the 
cut Ax^/(Xpeak/^dof) > 350. As a result, set A has 
fewer events removed by the supernova cut. 

In Alcock et al. (2000a) we describe in detail a 
robust way of statistically correcting for the t bias 
induced by blending. Briefly, this method is an inte- 
gration over the LF of the median t bias induced by 
blending in our sample of Monte Carlo events. As a 
check that this method gives a truly unbiased optical 
depth estimate we ran a series of secondary Monte 
Carlo simulations that make use of this correction 
and a number of Galactic halo models. The statisti- 
cal correction, although it blurs the individual events 
together, does a satisfactory job of reproducing an 
unbiased optical depth. 

The primary shortcoming of the present efficiency 
analysis is that all simulated events are assumed to 
be "normal" microlensing events with a single lens, a 
point source, and constant velocities. This assump- 
tion is used in the present analysis for simplicity and 
because of the highly uncertain characterization of 
exotic lensing events. A careful study using selection 
criteria much looser than criteria A & B has convinced 
us that it is unlikely that we have missed any exotic 
lensing events in the present data set. The primary 
concern is what effect the addition of exotic lensing 
might have on our detection efficiency (in particular 
because of binary lensing events). Although this is as 
yet uncertain it is probably a small effect due to the 
small number of exotic lensing events seen so far and 
the fact that criteria B does find some exotic lensing 
events. 

5. Event Distributions 

There are a number of statistical tests that can 
be performed on microlensing event distributions to 
test the hypothesis that events are gravitational mi- 
crolensing, or to test hypotheses regarding the lens 
population. As the sample of events becomes larger 
these tests become important discrimination tools. 

5.1. Impact Parameters 

An important model-independent test of the hy- 
pothesis that we have observed gravitational microlens- 
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ing is to compare the distribution of peak magnifica- 
tions to the theoretical prediction. It is convenient 
to switch variables from the maximum magnification 
(^max) to the minimum distance of approach between 
the MACHO and the line of sight, in units of the 
Einstein radius, Umin — b/rE- Events should be uni- 
formly distributed in Umm] this distribution is then 
modified by the experimental detection probability 
which is typically higher for small Umin ( Alcock et al. 



20C Oa ) . The observed and predicted distributions for 
our LMC events for both selection criteria are shown 
in Figure |[ A K-S test shows a probability of 79.8% 
of getting a K-S deviation worse than the observed 
value 0.172 for criteria A by chance, and a probabil- 
ity of 48.6% of being worse than 0.202 for criteria B. 
The binary event 9 is excluded from this comparison. 
We conclude that the distribution of events in Mmin is 
consistent with the microlensing interpretation. 

The Umin distribution and the high magnification 
events may be used to lend support to the microlens- 
ing interpretation of our lower magnification events. 
Our high magnification events are striking, and are 
clearly separated from the background in Figure |^. 
If these high magnification events are accepted as 
microlensing then there must exist many more mi- 
crolensing events with smaller peak magnification. 
Figure ^ shows that we find these in just the right 
proportion. 

5.2. The Color-Magnitude Distribution 

The gravitational microlens does not distinguish 
between types of stars, so naively one expects mi- 
crolensing to occur uniformly on every type of source 
star. However, both selection criteria A & B em- 
ploy various signal-to-noise cuts that bias us against 
detecting microlensing events on faint stars. In ad- 
dition, the measured baseline fiux of an event may 
be significantly larger because of the blending of non- 
lensed flux, and it is not always possible to accurately 
determine the amount of blending. That is, the blend- 
ing fits in Table ^ may not be reliable, since there is 
considerable fit parameter degeneracy between Amax 
and the blend fraction. Thus, detected microlensing 
candidates, while occurring on many types of source 
stars, may not follow the observed distribution of 
stars (or rather objects) exactly. 

Figure |^ shows a CMD with each of the 17 mi- 
crolensing candidates, along with all the closest 200 
stars around each candidate. Most of the events lie 
along the faint main sequence where most of the ob- 



served LMC stars reside. Event 5 is quite red for its 
brightness and could represent a foreground popula- 
tion of bright M-dwarf lenses as noted in A97, but as 
a whole the distribution in the CMD is quite repre- 
sentative. The distribution of events is not clustered 
in luminosity or the color magnitude diagram and is 
consistent with the microlensing interpretation. 

5.3. Spatial Distribution 

For microlensing by MACHOs smoothly distributed 
in the Galactic halo, or stellar lensing by stars in the 
Milky way thin or thick disk, we expect the detected 
events to be distributed across our fields in proportion 
to the local exposure. An extended LMC halo popu- 
lation could also form a smooth distribution. In con- 
trast, models in which LMC disk and bar stars dom- 
inate the lensing population predict that the lensing 
events will be concentrated within the LMC ((Alcock 
et al. 1997a; Aubourg et al. 1999; Salati, et al. 
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Gyuk et al. 1999), but also see Alves & Nelson 1999 



where a flared disk could widen this distribution to 
some extent). 

Figure |l| indicates that the detected events are ap- 
parently spread evenly across our 30 fields. To quan- 
tify this impression we perform two simple tests. For 
criteria A & B we computed a concentration param- 
eter, 9, as described in Gyuk et al. (1999). This 
parameter is a mean spatial separation between all 
combinations of events. For our two selection criteria 
we find 0A = 1.94° ± 0.29 and 9b = 1.86° ± 0.23, 
where the error bars have been estimated using the 
observed number of events and the models of Gyuk 
et al. (1999) (private communication). These num- 
bers should be compared with predictions from the 
various models of LMC self-lensing. Gyuk et al. 
(1999) find 9 = 1.3° ± 0.2 over the MACHO 30 fields 
for all their LMC disk-|-bar self-lensing models, and 
9 = 1.85°±0.15 for LMC disk-hbar-hhalo and Galactic 
halo models. Thus by this measure our event distri- 
bution is inconsistent with their most favored LMC 
disk-|-bar self-lensing at the ~ 2 cr level, but is con- 
sistent with an extended lens population such as is 
expected for a Galactic or LMC halo population. 

Our second test is illustrated in Figure ^ where 
the cumulative distribution in spatial distances on 
the sky, as measured from the optical center of the 
bar (a = 5/i24TO, (S = -69°48'), is plotted. Also shown 
are two predictions based on the models of Gyuk, et 
al. (1999): the dashed line is the predicted distribu- 
tion for uniform lensing (LMC disk-l-bar+halo) over 
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the face of the LMC, convolved with our detection 
efficiency per field (Alcock et al. 2000a), and the 



their masses or velocities. Experimentally, one can 
obtain an estimate of the optical depth as 



dotted line is the favored LMC disk+bar self-lensing 
model, also convolved with our detection efficiency 
per field. For criteria A, a K-S test yields a probabil- 
ity of 7.3% of getting a K-S deviation worse than the 
observed value 0.342 for the disk-|-bar model. For the 
disk-|-bar-|-halo model there is a probability of 59.5% 
of getting a K-S deviation worse than the observed 
value 0.204. For criteria B, a K-S test yields a prob- 
ability of 2.4% of getting a K-S deviation worse than 
the observed value 0.349 for the disk-|-bar model and 
for the disk-|-bar-|-halo model there is a probability of 
34.1% of getting a K-S deviation worse than the ob- 
served value 0.220. We note that these results are de- 
pendent on a direct comparison to the 30 fields used in 
this analysis, and that the LMC disk-|-bar self-lensing 
predictions are inconsistent with the data at the 93% 
C.L. 

One difference between a relatively smooth distri- 
bution of events over the face of the LMC caused by a 
Galactic population and that caused by a LMC halo 
is a slight east-west asymmetry due to the ~ 30° 
tilt of LMC's disk. This tilt would induce a slightly 
higher optical depth on the western side of the LMC, 
due to the longer path length through any LMC halo. 
No such east-west asymmetry should exist if the mi- 
crolenses are due to a Galactic population. A casual 
inspection of Figure Q hints at the possibility of such 
an east-west asymmetry. However, a simple Monte 
Carlo shows that the asymmetry is not statistically 
significant given the number of events in our sample. 

6. Implications 

We start with the implied microlensing optical 
depth, which is compared with the optical depth ex- 
pected from known populations of stars along the line 
of sight to the LMC. We then discuss our likelihood 
estimate of microlensing rate, MACHO masses, and 
optical depth for both the dark halo and known stellar 
populations. 

6.1. Optical Depth Estimates 

The simplest measurable quantity in a gravita- 
tional microlensing experiment is the microlensing op- 
tical depth, r, which is defined to be the instanta- 
neous probability that a random star is magnified by 
a lens by more than a factor of 1.34. This probability 
depends only on the density profile of lenses, not on 
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where i? = 6.12 x lO'' object-years is the total expo- 
sure, t i is the Einstein ring diameter crossing time 
of the ith event, and £{ti) is its detection efficiency. 
Here, and below, we use the statistically corrected val- 
ues of the blended durations t (Table |^). These take 
into account the fact that our typical star is blended, 
and so the fit t is typically underestimated. This sta- 
tistical correction depends upon the selection criteria 
used and is described more fully in § ^. It is also 
convenient to define the function 
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which is the contribution to Tmoas from a single ob- 
served event with timescale t . ti values for each of 
our events are also listed in Table |^ for both selection 
criteria. 

Using the criteria A set of 13 events, we find (Ta- 
bles H and ^) an optical depth for events with dura- 
tions 2 days < ? < 400 days of t|°° = 1.1 x 10"'^. 
With the criteria B set of 17 events, we find t^^" = 
1.3 ^p'g X 10~^. This is to be compared with r = 
4.7 X 10~^ for a typical dark halo consisting entirely 
of MACHOs and with predicted Tatars = 0.24 x 10""^ 
to 0.36 X 10~^ from known stellar populations (from 
Table ^ below) . Subtracting the stellar lensing back- 
ground from our observed optical depth, we find that 
the observed excess is about 15%-25% of the predicted 
microlensing optical depth for a typical all-MACHO 
halo of equation (jj) below. 

This optical depth estimate has the virtue of sim- 
plicity; however, since the events are "weighted" oc ri, 
it is hard to assign meaningful confidence intervals to 
T without assuming some particular t distribution. 
This is illustrated in Figure |[ which shows the contri- 
bution to the sum in equation (|l]) from events in var- 
ious bins oft . In comparison with A97, we note that 
while the contributions become large at large t , they 
are substantially smaller in the 100-300 day range due 
to the increased baseline and looser cut t < 600 days. 
The large contribution at long t implies that the over- 
all uncertainty in r is greater than simple Poisson 
statistics based on 13 or 17 events. However, this un- 



14 



certainty will continue to decrease as the experiment 
progresses. For example, in A97 we estimated that 
should we expect to have observed on average 1 ad- 
ditional event with t ^ 300 days (but happened to 
observe no such event), the real r would have been a 
factor of two higher, and we were not able to exclude 
such a possibility with any confidence. The equivalent 
situation with the current data set is less dangerous, 
as such a missed event would result in a real t being 
only about 20% higher. However, our optical depth 
estimate is valid only for a specific mas s or timescale 
interval. The likelihood method of § 12 gives another 
way of estimating the optical depth and the errors on 
the optical depth. 

6.1.1. Optical Depth Cut Dependence 

Figures |l^ and |l^ show the dependence of the the 
measured optical depth on the Umin and Ax^/ (Xmi/^dof ) 
cuts. The heavy curves indicate Tmeas for set A while 
the light curves give Tmeas for set B. For the binary 
event we have assigned a Umin value of 0.573, which 
is the value obtained for the single-lens fit. 

Figures |l^ and ^ clearly indicate that our opti- 
cal depth results are not very sensitive to our cut 
values. The Tmeas values generally do not vary by 
more then the Icr statistical error bars for Umin and 
Ax^/(Xmi/^dof) cuts in the ranges 0.2 < Wmin < 1-0 
and 300 < AxV(Xmi/^dof) < 5000. We note the 
largest single contribution is from event 22, which was 
included in set B and excluded in set A. 

6.1.2. Comparison of A97 with the Present Analysis 

Why is our new value of the optical depth a factor 
of two smaller than the value reported in A97? The 
reasons are manyfold and somewhat difficult to sepa- 
rate out completely. By far the largest effect can be 
classified as Poisson in nature. We list the causes in 
order of decreasing effect on the optical depth. 

Inspection of Figure |^ reveals that a disproportion- 
ate number of our 'high' quality events were observed 
in the first 2.1 years of the data set. Events 1, 5, 7, 
and 9 are all of high quality. Only event 14, 21, 23, 
in the following 3.6 years stand out as having compa- 
rable quality. This 'qualitative' feel for the events is 
backed-up by the fact that in the first 2.1 years, 6-7 
events were observed: a rate of 2.9-3.3 events/year. 
And for the last 3.6 years only 7-10 events were ob- 
served: a rate of 1.9-2.8 events/year. Since the effi- 
ciency has not changed drastically over this interval. 



we conclude that while the exposure increased by a 
factor of 3.4, the number of events did not. That 
is, we got "lucky" during the 2.1 years of A97 and 
detected more microlensing than average. This possi- 
bility was reflected in the large Poisson errors quoted 
in A97. 

Another reason the optical depth reported in A97 
is large compared to this work is related to the split- 
ting of the top 6 fields for analysis described in § |[ 
These 6 fields represented ^ 27% of our exposure in 
object-years. Because they were split into two sets 
of lightcurves of approximately one year duration the 
efficiency was lowered for durations > 100 days. At 
least 4 of the events in the 8 event sample of A97 
(events 5, 6, 7, 9) contributed somewhat more to the 
optical depth then they would have had the top 6 
fields not been split for analysis. This in itself was 
not in error, for one would have expected that one 
or two events would have been missed because of the 
splitting of the top 6 fields and these 'missed' events 
would have counter-balanced the decreased efficiency. 
By chance, no such events were missed, as is evident in 
the current analysis where the top 6 fields have been 
merged and analyzed as full lightcurves. We estimate 
that this "small numbers" (Poisson) effect increased 
the optical depth reported in A97 by ~ 25%. 

For reasons discussed in more detail in § ^ (and 
Alcock, et al. 2000a) , our detection efficiency is some- 
what higher than in A97. Briefly, we neglected a con- 
tribution from faint stars with V > 21.5 and thus 
our efficiency in A97 had not yet converged for dura- 
tions > 100 days. In addition, our new normalization, 
which has been more carefully determined using HST 
data, and has more realistically estimated errors leads 
to a somewhat increased sensitivity. Together, these 
effects spuriously increased the optical depth in A97 
by ~ 10% with respect to the current results. This is 
within our estimated uncertainty. 

As mentioned in § 3.3.2, one of the events used in 
A97 (event 10) is most likely a SN interloper. This 
interpretation is supported by the presence of an ob- 
vious spiral galaxy in our HST frames of this event 
and the fact that it is quite reasonably fit by a type 
la SN lightcurve. The effect of this interloper spuri- 
ously increased the optical depth by ~ 7% in the 8 
event sample and had no effect on the 6 event sample 
(which rejected it). 
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6.1.3. Errors in Present Analysis 

While the current analysis is the most careful yet 
performed and a substantial improvement over earlier 
efforts, there are still a number of errors or potential 
errors that exist in our results. The errors due to 
small number statistics are included in the error bars 
we report (about 30% uncertainty). The errors due 
to model dependency are explored by considering a 
range of models. We believe these are the largest 
errors in our results. In this short section and in Ta- 
ble |l^ we list some other sources of error that might 
be worth considering in more detail in the future. 

Normalization of star to object ratio: The ratio of 
actual LMC stars to SodoPhot objects varies across 
our fields, and it is very difficult to accurately esti- 
mate. We have HST images for three areas, and at- 
tempt to tie together MACHO object based photom- 
etry to the HST star-based photometry to create a 
unified luminosity function (LF). We estimate a 20% 
uncertainty in our final results because of this. Is- 
sues include the underlying luminosity function, the 
magnitude calibration of our objects, blending effects 
in matching the object LF to the stellar LF, and the 
unknown effects of crowding, seeing, and sky in the 
template images, among several others. 

Selection Criteria: Since we do not have a com- 
plete understanding of the background, and because 
we examine our data before deciding upon the selec- 
tion criteria, it is possible that we differentiate back- 
ground and microlensing in an imperfect and/or bi- 
ased way. We implemented two independent sets of 
selection criteria as a test of our sensitivity to this 
bias. We estimate about a 20% uncertainty due to 
our selection criteria. 

Correction to t : Blending causes the durations of 
events to shift from their naive fit values. We chose to 
correct for this in a statistical manner, and estimate 
about a 3% uncertainty in this correction. This is 
due to several factors, but mainly the uncertainty in 
the true t distribution which is needed to make the 
correction. 

Binary source stars: Locally most stars reside in 
binary or multiple star systems, and it is expected 
that this is also true of LMC stars. Our LMC lumi- 
nosity function does not include a correction for this. 
This is a complicated correction which will be uncer- 
tain since the binary fraction for LMC stars is not 
known. We did not make an estimate of the size of 
this effect, but will consider it in a later paper. 



Exotic microlensing: We do not explicitly add 
binary-lens microlensing or other exotic lensing into 
our artifical lightcurvcs. Thus our efficiencies for 
these are unknown. We did try to explore two sets of 
selection criteria to detect any gross sensitivity of our 
results to this effect, but a proper calculation should 
be done in the future. 

Others: There are several other sources of system- 
atic error that have been considered. For example, 
in the artifical star Monte Carlo we assume that all 
fiux added goes to the nearest SoDoPhot object. Di- 
rect tests show that this is not true in about 3% of 
the cases. The effect of this mis-identification is not 
known, but could be around 3%. Some other small 
errors are discussed in Alcock, et al. (2000a) and Van- 
dehei (2000), and a summary is given in Table |l^. 

Finally we note that due to the complex nature of 
systematic errors, it is not straightforward to combine 
them with each other or with statistical errors. For 
example, systematic errors can be strongly asymmet- 
ric. 

We estimate our total systematic error to be in the 
range of 20% -30%, though even this range is uncer- 
tain. 

6.2. Likelihood Analysis and Dark Matter 

We compare the number of detected events and 
the distribution of observed timescales, t, with pre- 
dictions from models of various lens populations. Mi- 
crolensing can occur when any compact object travels 
in front of a monitored star, so we expect microlens- 
ing events from any population of stars, remnants, 
or dark compact objects that lie between us and the 
LMC. Luckily, much is known about the density and 
velocity distribution of stars and remnants in the 
Milky Way and LMC. Less is known about the dark 
halo of the Milky Way (and even less about the dark 
halo of the LMC), but we can leave the fraction, /, of 
dark objects that can lens, as well as the masses, m of 
these objects, as free parameters which we determine 
using a maximum likelihood analysis. This analysis 
was done in A97, with 6 or 8 observed microlensing 
events. We found values of / between 0.15 and 1.0 
at 90% confidence level (CL), and MACHO masses 
m between 0.1 and 1.0 Mq. The large uncertainty in 
these results came largely from small number statis- 
tics, but also from uncertainty in the models. In the 
current analysis, we have a larger number of events 
and an improved efficiency determination, so we can 
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reduce the Poisson and some systematic errors. Here 
we also make improvements to our likelihood analy- 
sis, most notably the inclusion of realistic estimates 
of the stellar lensing background within the likelihood 
function. We note that the results of this section still 
depend heavily upon the models of Milky Way and 
LMC. We will come back to this point in § |^. 

6.2.1. Microlensing Rate 

The microlensing event rate T is more model- 
dependent than the optical depth r. Rate F depends 
on the event timescales via the mass function of MA- 
CHOs and their velocity distribution, but the uncer- 
tainties in r are given purely by Poisson statistics. 
Thus F is useful in quantifying the errors on any mea- 
surement, once a halo model is specified. 

The number of observed events is given by a Pois- 
son distribution with a mean of 
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For a typical dark matter halo (equation consist- 
ing of 100% MACHOs, the total rate of microlensing 
events with Amax > 1-34 is given by equation (A2) 
of A96: F w 1.6 X W-^{m/ Mq)-°-^ events/star/yr. 
Thus if all MACHOs had the same mass and our effi- 
ciency were 100%, we would expect about 100(to/ Mo)^^-^ 
events in the present data set. The average timescale 
of an event scales oppositely, (i) ~ 140(m/ Mq)^/^, 
since the product of the two gives the optical depth 
which is independent of the MACHO masses, r = 
■|F(i). Thus, although the optical depth is inde- 
pendent of MACHO mass, in a real microlensing ex- 
periment statements on the MACHO content of the 
halo will be dependent on the MACHO mass. How- 
ever, the masses of the lenses are constrained, since 
we measure the distribution of event timescales. 

6.2.2. Milky Way and LMC Models 

We consider four stellar components and a dark 
component of the Milky Way, and a stellar and dark 



component of the LMC. Given our exposure and effi- 
ciency, and a model of the density, velocity distribu- 
tion, and mass function of a lens population, we can 
calculate the expected microlensing optical depth, mi- 
crolensing rate, distribution of event durations, and 
the number of expected events detected in our exper- 
iment. A summary of the results are shown in Ta- 
bles |l^ and |l^. For the LMC self-lensing model, we 
note that the rate depends strongly on the position 
on the sky, so the values we report depend on the 30 
specific fields we monitor. 

We model the density of the Milky Way and LMC 
thin and thick disks as double exponentials 
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where z and R are cylindrical coordinates, Mdisi? is the 
total mass of the disk, Zd is the scale height, and Rd is 
the scale length. Instead of specifying the total mass. 
Milky Way disks are often specified by the column 
density, Eq, at the solar circle, R = Rq ^ 8.5 kpc, 
and the relation is 



2nRi 

For simplicity, we characterize the velocity distribu- 
tion of a disk as a constant rotation velocity, Vc, with 
some isotropic dispersion, a^, in addition. 

For the normal Milky Way thin disk we use param- 
eters: Rd — 4: kpc, Zd — 0.3 kpc, Sq = SOMqpc"'^, 
Rq = 8.5 kpc, Vc = 220 km/s, and cr^, = 31 km/s. 
This gives a total thin disk mass of Mdisk = 4.2 x 
10^" Mq. Later, we will discuss models (e.g., model 
F) with a maximal thin disk, and a smaller dark halo. 
In that case we use the above parameters except with 
So = 8OM0pc-3^ which gives Mdisk = 6.7 x 101°Mq. 

For the Milky Way thick disk we use parameters: 
i?d = 4 kpc, Zd = 1.0 kpc. So = 4Mqpc-3, Rq = 8.5 
kpc, Vc = 220 km/s, and — 49 km/s, for a total 
mass of Mdisk = 3.4 x 10^ M©. 

For the LMC disk we use the preferred param- 
eters from Gyuk et al. (1999), Rd — 1.57 kpc, 
Zd = 0.3 kpc, Vc = 70 km/s, — 25 km/s, and 
Mdisk = 3.0 X IO^Mq. The LMC disk self-lensing 
also depends upon its distance, L = 50 kpc, inclina- 
tion, i — 30°, and position angle, (f) — 170°. These 
are the parameters we use when considering the LMC 
disk plus LMC dark halo model. For the pure disk (no 
LMC halo) case, we conservatively increase Mdisk to 
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5 X 10^ Mq, in good agreement with recent anal ysis of 
the LMC rota tion curve (Mdisk = 5.3± I.OMq ( |Alves 
et il. 1999d), corresponding to central surface den- 
sities of 190 Mqpc"^ and 320 M0pc~^ respectively. 
We do not consider a separate bar component, since 
the bar mass is strongly limited by the HI kinematics 
dKim et al. 1998|; pyuk et al. 1999)). 



The Milky Way spheroid density is modeled as 
dGuidice et al. 1994|; [Alcock et al. 1996aD 
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with no rotation, and an isotropic velocity dispersion 
of cTt, = 120 km/s. 

For the LMC we consider two main cases: 1) pure 
disk/disk self-lensing, and 2) disk/disk self-lensing 
plus an LMC dark halo consisting of fraction / of 
MACHOs, with / being the same fraction used for 
the Milky Way dark halo. Later, we also consider the 
possibility of an all stellar LMC halo. 

For the Milky Way dark halo we consider three 
models: S, B, F, which were used in A96 and A97. 
The density of Model S is given by 



PH{r) = po 
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where pn is the halo density, po — 0.0079 Mq pc~^ 
is the local dark matter density, r is Galactocentric 
radius, i?o = 8.5 kpc is the Galactocentric radius of 
the Sun, and a — b kpc is the halo core radius 
With the standard thin disk, this model has a to- 
tal rotation speed at 50 kpc of 200 km/s, with 190 
km/s coming from the halo. We assume an isotropic 
Maxwellian distribution of velocities with a 1-D rms 
velocity of 155 km/s, and assume a (5-function MA- 
CHO mass function of arbitrary mass m. ^ Note 
that we always multiply the above density by the MA- 
CHO halo fraction /, implicitly assuming that the re- 
maining 1 — / fraction of the halo is filled with some 
exotic particle dark matter or other non-lcnsing mat- 
ter. Dark halo models B and F are power-law models 
(Evans 1993; Alcock et al. 1995d), and are discussed 



in detail in A96, and A97. Model B has a very large 



^'^ Analysis of carbon star kinematics on the periphery of the LMC 
disk support a, psniido-isothermal density profile for the Galac- 



tic dark halo ( Alves et al. 1999d| ). 

As discussed in A97, a delta-function mass distribution is a rea- 
sonable fit to the observed distribution, and more complicated 
forms are difficult to distinguish with such a small number of 
events. 



dark halo and a standard thin disk, giving a rising ro- 
tation curve that reaches 258 km/s at 50 kpc. Model 
F has an extremely low mass halo, somewhat incon- 
sistent with the known Galactic rotation curve, and 
a very massive disk. At 50 kpc, the model F halo 
contributes 134 km/s towards a total of 160 km/s ro- 
tation speed. 

Finally, for the LMC halo we follow Gyuk et al. 
(1999) and use the same density distribution as model 
S above, but with a central density of 0.0223 Mqpc"'^, 
a = 2 kpc, Vc = 70 km/s, and a tidal truncation 
radius at 11 kpc. In this model the mass of the LMC 
halo in the inner 8 kpc is 6 x 10^ M© and the total 
mass of the halo is 9.2 x 10^ M©. This is a somewhat 
extreme model, probably larger than allowed by the 
LMC rotation curve. Like the Milky Way halo, the 
LMC halo is assumed to consist of a fraction / of 
MACHOs all of mass to. It is conceivable that the 
LMC MACHO halo fraction and make-up is different 
than the MACHO fraction in the Milky Way halo, 
and we consider this possibility in the next section. 
We note that no substantial stellar component of an 
LMC halo has yet been observed. 

For the stellar lensing populations we integrate the 
microlensing rate over a mass function. There have 
been several recent determinations of the present day 
mass function (PDMF), but it is not clear that the 
mass function determined locally is valid for all the 
stellar populations we model. However, for simplicity 
we will use the PDMF of Gould, BahcaU, & Flynn 
(1997) for all the stellar populations. Table |ll| shows 
some properties of the stellar population calculated 
from the models above. It also shows the expected 
number of microlensing events from each population 
found by the likelihood method, and thus including 
the effect of our efficiency calculation and selection 
criteria. These results differ to some degree from 
those presented in A97 for several reasons. First, the 
models we use are different in some cases, and we are 
using a different PDMF. Second, and most impor- 
tantly, we explicitly do not count lenses that are too 
bright to be detected as microlensing in our experi- 
ment. We have an explicit cut at around V = 17.5 
magnitudes, and so stars brighter than this cannot 
be found as lenses. Our Monte Carlo shows that this 
cuts the expected number of thin disk stellar lensing 
events by more than half, with smaller effects for thick 
disk, spheroid, etc. and almost no effect for LMC disk 
lensing. The results displayed in Table |ll|, use the full 
Monte Carlo for all stellar distributions. 
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6.2.3. MACHO Halo Fraction and Mass 

We find the most likely values of halo fraction / 
and unique MACHO mass m using our set A (13 
events) and set B (17 events), and their corresponding 
efficiencies. The likelihood function is 



L(m, f) = exp(-iVcxp) 
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theoretical rate of microlensing derived from model 
j. The distributions ^(«) for stellar populations are 
integrated o ver a mass function ( [Gould, Bahcall,"fc 
Flynn 1997 ) as described above, and are calculated 
using code described in Gyuk et al. (1999). The 
results are dependent on the model so we explore 
a range of possible halos, including a standard halo 
(model S from A96 and A97). We also use two other 
dark halo models. We choose model B from A96 and 
A97, because it is about as large a halo as the data 
will allow. We also choose model F from A96 and 
A97 because it has a nearly maximal disk and a very 
low mass halo, and therefore is as small a halo as the 
data allows. These models are described in detail in 
A96. We do not show models A, C, D, or G from A96 
and A97, or other possible halo models since they are 
in general intermediate between the extremes of mod- 
els B a nd F. Model S is a common pseudo-isothermal 
sphere ( Griest 1991 ) with an asymptotic rotation ve- 
locity of 220 km/s, while model s B and F are power- 
law Evans models ( Evans 1993 ). Tabl e [l^ sh ows the 
results for all the models, and Figures |l2|, |13|, and |lj 
show the corresponding likelihood contours. 

For model S, the resulting likelihood contours, as- 
suming a (5-function halo mass function, are shown 
in Figure 12, The probabilities are computed using a 
Bayesian method with a prior uniform in / and log m. 
We show likelihood contours for both our 13 event 
sample and our 17 event sample, and with and with- 
out the LMC halo. The best fit values and errors are 
given in Table |lj. The errors shown in the table are 
one-sigma errors. The peak of the likelihood contours 
gives the most probable mass and halo fraction for a 
given model and for set A with an LMC halo we find 
TO2D = 0.48 Mq, /2d = 0.20. For the corresponding 
set B, m2D = 0.67 M©, and /2D = 0.23. 



We calculate the one-dimensional likelihood func- 
tion by integrating over the other parameter and 
find (for set A without an LMC dark halo) a most 



likely MACHO mass toml = 0.60^0 20 Mq, and most 
likely halo fraction /ml = 0.211°;^°. The errors 
given are 68% CL. The values for set B are tuml = 
0.79t°J^MQ and /ml = 0.24^° ;;^. For model S, the 
95% CL contour includes halo fractions from about 
8% to about 50%, and MACHOs masses from about 
0.12 Mq to 1.1 M0, depending upon the selection cri- 
teria and LMC model used. The likelihood method 
gives an optical depth for the halo population of 
l.lj^Q'4 X 10^^ almost independent of the section cri- 
teria, the LMC model, and the Galactic model. 

There are several important comments to be made. 
First, sets A and B give results that are remark- 
ably similar, implying that the systematic error intro- 
duced by our selection criteria methodology is small. 
The important parameters of estimated MACHO halo 
fraction are nearly identical using the two different 
sets of events and efficiency determinations. The es- 
timated typical MACHO mass does vary between the 
two sets of events, but the values lie within one sigma 
of each other. This difference in lens mass comes par- 
tially from the rejection of event 22 from set A. 

Second, consistent with our optical depth esti- 
mates, the values of the halo fraction are approxi- 
mately a factor of two lower than we found in A97. 
As discussed in § 6.1.2 , this is mainly a result of find- 
ing more events per unit exposure during the first 
two years, but is also due to changes in efficiency, etc 
We note that the optical depths reported in Table 12 
are the estimated MACHO contribution, and do not 
include the background of stellar microlensing. The 
contributions from stellar background are shown in 
Table ^ The values found here are quite similar to 
those found directly in § 3.1 
Third 



our new confidence intervals are substan- 
tially smaller than those of A97 due to the larger 
number of events. Even though the central values 
have changed, there is reasonable overlap of our new 
contours with the A97 contours, and our new most 
likely values lie within the A97 90% confidence re- 
gion. 

Next, for Model S with a large LMC disk, but no 
LMC dark halo, and set A, we expect a total of 3.0 
events from stellar background sources, with the ma- 
jority coming from LMC self-lensing. For the same 
model and set B, the number of expected background 
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events is 3.9. In both cases the predicted number of 
background events is substantially below the number 
of detected events. Thus, if these models are correct, 
the microlensing events are very unlikely to come from 
the known stellar populations. 

For the case of an LMC halo plus LMC disk, LMC 
disk self-lensing must be smaller since part of the 
LMC rotation curve is supported by the halo. In 
this case some of the lensing can come from the dark 
halo. This changes the predictions of MACHO halo 
fraction, since the LMC halo contributes very little 
to the total mass of the Milky Way, but relatively 
more to the microlensing. As shown in Table 
for Model S we find 1.1 events from the LMC halo, 
and 2.1 background events using the set A. For set 
B we find 1.4 LMC halo events, with 2.7 background 
stellar events. Again, the expected number of back- 
ground events is significantly smaller than the number 
of observed events. When an LMC dark halo is in- 
cluded the events from the LMC halo count towards 
dark matter that is not uniformly spread across the 
sky. The predicted values change from / — 0.21 to 
/ = 0.20 and from m = 0.60 Mq to m = 0.54 Mq 
for set A, and from / = 0.24 to / = 0.22, and from 
m = 0.79 Mq to m = 0.72 M© for set B. The change in 
MACHO fraction is small because our LMC halo has 
an optical depth of 0.79 x 10~^, substantially smaller 
than the 4.7 x 10^^ contributed by Model S, and the 
LMC disk contributes more background when no halo 
is present. The most likely total mass in MACHOs 
in the Milky Way dark halo (within 50 kpc) shows an 
expected drop of about 10%, (from 8.5 x 10^" M© to 
7.9 X IO^^Mq) when an LMC dark halo is included. 
As we discuss below, with Model F, which has a very 
small Milky Way halo, the change in / and the change 
in total MACHO mass is much more substantial. 

Finally we note that with a typical halo model like 
S, the likelihood contours in all cases rule out a 100% 
MACHO halo at high significance. This was not true 
in A97, or in any previous microlensi ng data set. (B ut 
see the very recent EROS II reports ( Milsztajn 1999 )). 
We note that our sensitivity to events longer than a 
thousand days is small, so we cannot rule out dark 
matter objects with masses of tens of solar masses. 

6.3. Interpretation 

Several interesting features can be seen in Table [ij 
and Figures |l|, [l|, and |lj. 

Examination of the likelihood contours show that 



with our new data set, the uncertain nature of the 
Milky Halo dominates over Poisson error and the sys- 
tematic error caused by our selection criteria. For 
each model, the most likely values and confidence lim- 
its are nearly the same, while between models there 
are significant differences. This is an improvement 
over A96 and A97 where small number statistics dom- 
inated the errors. The values found are typically 
within 1.5 sigma of those reported in A97, but a factor 
of two smaller for reasons given above. 

As noted in A96 and A97, the most likely halo 
fraction, /, depends strongly on the halo model, with 
massive halos such as Model B giving small MACHO 
fraction (/ ~ 13%), medium halos such as Model S 
giving medium values (/ ^ 22%), and very low-mass 
halos such as Model F giving large fractions (/ ^ 40% 
- 60%). However, there are some model independent 
conclusions that can be drawn. The total predicted 
mass in MACHOs within 50 kpc (column 5) is about 
9 X 10^*^ Mq for all models. This is again a factor of 
two smaller than reported in A97 for the same rea- 
sons. However, for very small halos, such as model 
F, the total MACHO mass is somewhat dependent 
on the model of the LMC halo. With no dark LMC 
halo, masses up to 10 x 10^° M© are found, while 
with a large dark LMC halo, the prediction drops to 
8 X 10^° Mq. This is because this MW halo has an op- 
tical depth of only 1.9 x 10^ compared to 0.79 x 10^^ 
for our LMC halo. So with the large MACHO fraction 
caused by the very small MW halo, the LMC halo can 
contribute substantially to microlensing without con- 
tributing much to the mass within 50 kpc. Note, how- 
ever, that when one sums the total MACHO contri- 
bution to microlensing optical depth (MW halo plus 
LMC halo), the result is tml = 1-1 ± 0.4 x 10""^ al- 
most completely independent of the MW and LMC 
halo models. This is shown in column 6 of Table |l2|. 
These values and their confidence intervals are simple 
to interpret statistically, since each model provides a 
distribution of event durations. Thus the subtleties 
discussed in Section 6.1 are absent. We note that 



these values include only the LMC and MW halo con- 
tributions, and are close to the values we obtained in 
our direct estimates of optical depth. 

Interestingly, for models S and B, halos consist- 
ing 100% of MACHOs are strongly ruled-out. Even 
for the rather extreme model F, and no LMC halo, 
a 100% MACHO halo is ruled-out. The only way 
this data is consistent with a 100% MACHO halo is 
if there is an extremely small MW halo coupled with 
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a very small LMC halo. This is the strongest limit to 
date on an all MACHO halo and is a major result of 
this work. We note that we do not set strong limits on 
dark matter objects with masses in the tens of solar 
mass range. However, the 9 x 10^° Mq in MACHOs 
found in this work still represents several times the 
mass of all known stellar components of the Milky 
Way. If the bulk of the lenses are located in the halo, 
then they represent the dominant identified compo- 
nent of our Galaxy, and a major portion of the dark 
matter. 

While tmacho f^nd the total mass are fairly model 
independent, as discussed in A96 and A97, the typi- 
cal mass of a MACHO, m, is not. Lighter halo mod- 
els such as F have a smaller implied MACHO mass 
m-ML ~ 0.2 Mq, while heavier halos such as B have 
mML 0.8 Mq, and medium halos such as S giving 
m-ML ^ 0.6 Mq. Sets A and B also differ slightly in 
this parameter. We conclude that our estimate of m 
is not very robust, but that masses below the brown 
dwarf limit of 0.08 Mq are quite unlikely. Therefore, 
the nature of the lenses remains unclear. 

Finally, using Table |l^, and comparing the num- 
ber of expected events from the MW halo with the 
number of expected events from known stellar pop- 
ulations, one can ask how strong is the case for any 
MACHO contribution to the dark matter. For set A, 
the known stellar background ranges from 2.1 to 3.2 
events depending upon the LMC and MW halo model. 
For set B, the range is 2.7 to 4.2 stellar events. If our 
models of the MW thin disk, thick disk, spheroid, and 
LMC disk are adequate then it is very unlikely to find 
13 (or 17) events when expecting no more than 3.2 (or 
4.2) (probability less than 10^^ in both cases). We 
note that the LMC disk we use in the pure-disk model 
is conservatively large (Mdisk = 5 x 10^ Mq). 

The question remains whether an LMC halo could 
supply the observed microlensing. The halo we use is 
a good fit to the available data, (Mhaio = 6 x 10^ Mq 
within 8 kpc), and the rotation curve for this halo 



plus LMC disk is a good fit to the data (Gyuk et al. 
19£|9|). Now, due to the lack of stellar tracers found 
with velocity dispersion of the ~ 50 km/s expected for 
a halo population, we have assumed that this LMC 
halo is dark. A large dark halo is, of course, typical of 
dwarf spiral galaxies such as the LMC. So assuming 
the halo consists of a fraction / of MACHOs, is a rea- 
sonable first approximation. Table |l2| shows almost 
no difference in predicted background or other quan- 
tities for Models S and B that have medium to large 



dark halos; however, for the very small halo Model 
F, nearly half the expected events come from stars or 
the LMC halo. The likely final halo fraction, is still 
nearly 50%, but with substantially larger error bars. 
Figure [l^ shows, however, that even in this 
no-MACHO halo is quite unlikely. 



Recently, several groups (Aubourg et al. 1999 
Weinberg 1999; Gyuk et al. 1999D have considered 



lensing by an extended stellar population around the 
LMC, and have obtained different estimates of the 
optical depth contribution. Interestingly, Graff, et al. 
(1999) claim tentative (^ 2<j) detection of a kine- 
matically distinct population. We therefore explore 
the effect this would have on our maximum likeli- 
hood analysis, by considering the effect of an LMC 
halo consisting of MACHO fraction / and a Milky 
Way halo with no MACHO population (thus presum- 
ably consisting entirely of exotic elementary parti- 
cles). The likelihood maximum for such a model gives 
/ = 1.35t5];t and m = 0.2^0:^8 M© for set A, and 
gives / = 1.52ti^;t and m = 0.31[!:o9Mq for set B. 
Thus a possible non-dark-matter explanation for our 
results is an LMC halo of mass ~ 9 x 10® Mq, consist- 
ing of stellar-like objects that have not yet been con- 
vincingly observed. This value of LMC halo mass is 
somewhat extreme, though there are published mod- 
els with masses this high. Clearly it is important to 
discover, or convincingly rule-out, the possibility of a 
large LMC stellar halo. 

7. Summary and Discussion 

We have detected between 13 and 17 microlensing 
events towards the LMC. The implied optical depth, 
microlensing rate, and MACHO halo fraction are a 
factor of two smaller than found in our previous work, 
but are consistent with previous results within the er- 
rors of small number statistics. The larger number of 
events allows us to reduce the Poisson error consider- 
ably, which along with our improved efficiency anal- 
ysis and study of sources of systematic error, means 
that the interpretation of the microlensing events is 
now dominated by uncertainties in the models of the 
Milky Way and LMC. We find that the number of 
events is not consistent with known lens sources and 
our measured optical depth, t|°'^ = l-2lo;3 x 10"'', 
is significantly larger than allowed by known Galactic 
and LMC stellar populations. The total implied mass 
in MACHOs within 50 kpc is ~ 9 x 10^° Mq, quite 
independent of the dark halo model. This is substan- 
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tially larger than all known stellar components of the 
Galaxy. However, one of our most important conclu- 
sions is that a 100% all MACHO Milky Way halo is 
ruled out at the 95% C.L. for a wide range of reason- 
able models. 

One explanation of our results is a Milky Way halo 
consisting of about 20% MACHOs. Another possibil- 
ity is an LMC halo that dominates the microlensing, 
and no MACHOs in the Milky Way halo. The spatial 
distribution of events makes lensing entirely by a stel- 
lar population in the LMC disk or bar less likely, but 
given the highly uncertain nature of an LMC halo, 
a previously unknown stellar component or an LMC 
halo consisting of MACHO dark matter, could explain 
the measured optical depth, the number of observed 
events and their spatial distribution on the sky. We 
note, however, that no known LMC stellar popula- 
tion exhibits kinematics of this nature. It would help 
to constrain the kinematics of old populations in the 
LMC and to look for new populations that could rep- 
resent the lenses. A direct measurement of the dis- 
tance to some LMC lenses would be especially useful 
in distinguishing the two possibilities above. 

There is intriguing evidence for a population of 
white dwarfs (Ibata et al. 1999; Mendez & Min- 
niti 1999) in the Hubble Deep Fields (North and 
South), consistent with the Milky Way halo hypoth- 
esis. These results were spawned primarily by the re- 
cent work of Hansen (1999) on low temperature white 
dwarf cooling curves. Although the identification of 
these faint-blue objects as white dwarfs remains to be 
confirmed and the small sample size restricts an accu- 
rate estimate, the suggestion that these white dwarfs 
could contribute 1/3 to 1/2 of the dark matter in 
the Milky Way is certainly stimulating in light of the 
present work. A third epoch HDF image to check the 
proper motion of these objects should help confirm or 
rule-out this hypothesis. A Galactic halo composed 
of ~ 20% by mass of white dwarfs would seem to be 
a natural explanation of both the microlensing data 
and this newly observed population, but the forma- 
tion of such objects and the chemical enrichment they 
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Table 1 
Field Centers 



J7 11, i( i _N n . 


V (. Ill <. J. . lA 1 \. 




V. J 1 I. \ l l\ iUlio 


1 


05 05 23 


-69 05 24 


1017 


2 


05 12 47 


-68 30 21 


860 


3 


05 22 24 


-68 28 01 


720 


5 


05 11 17 


-69 40 18 


839 


6 


05 20 00 


-70 17 10 


856 


7 


05 28 54 


-70 27 31 


1027 


9 


05 10 57 


-70 23 40 


811 


10 


05 04 34 


-69 52 19 


665 


11 


05 36 56 


-70 31 34 


930 


12 


05 45 36 


-70 35 16 


772 


13 


05 19 39 


-70 51 40 


752 


14 


05 35 53 


-71 09 22 


741 


15 


05 45 34 


-71 14 36 


718 


17 


04 57 04 


-69 43 11 


345 


18 


04 57 55 


-68 56 08 


594 


19 


05 06 09 


-68 21 03 


672 


22 


05 11 18 


-71 00 19 


398 


23 


05 02 51 


-70 35 39 


339 


24 


05 00 39 


-67 56 45 


297 


47 


04 53 05 


-68 01 26 


610 


53 


05 02 09 


-66 40 19 


180 


55 


05 02 15 


-65 58 50 


241 


57 


05 09 07 


-65 46 07 


189 


76 


05 44 13 


-69 49 41 


386 


77 


05 27 24 


-69 45 24 


1338 


78 


05 19 26 


-69 42 27 


1312 


79 


05 12 59 


-69 05 43 


1226 


80 


05 22 44 


-69 05 18 


1186 


81 


05 35 56 


-69 49 34 


792 


82 


05 32 50 


-69 03 18 


757 



This table lists the 30 well-sampled fields used in the 
current analysis, the 8 new fields arc 17, 22, 23, 24, 53, 55, 
57, & 76 and typically have less than half the observations of 
the top 22 fields. We observe 82 LMC fields in total, but the 
remaining 52 were observed less often (~ 120 observations 
each) . 
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Table 2 
Selection Criteria 



Description 



Year 2 (A97) 



Criteria A 



Criteria B 



Min. coverage 



bmrN >= 7 
> 40 baseline points 
t < 300 



SN87A echo 10' x 10' square excl. 

Crowd & PSF fcRD < (AxV(xV^dof)) /520 

& crdrej < 0.05 



Bumper cut 

Variable cut 
High points 
Baseline fit 

2nd S/N 
Main S/N 
Magnification 

2nd peak 
Supernova cut 



V > 17.5 
& y - < 0.9 



None 

6 pts. > 2(7 
& > 1 pt. on rise & fall 

Xml-out/^dof < 4 
AxV(X?eak/^dof) > 200 

AxV(xii/iVdof) > 500 
> max(1.75, 1 + 2a) 
None 
By eye 



riV > & &iV > 
> 45 baseline points 
t < 600 

^max ^ 310 

10' X 10' square excl. 
None 



V > 17.5 
& y - i? < 0.9 



bauto/rauto > 0.75 
7 pts. > 2a 

Xil-out/^dof < 1.8 
AxV(X?eak/^dof) > 350 

AxV(xii/iVdof) > 400 
>niax(1.49, 1 + 3.5ct) 
None 

AxsN-ML > & not event 22 



> 65 simul. 
baseline points 
t < 600 

^max ^ 310 

10' X 10' square excl. 

pkcrdrej + 
pkpsfrej < 0.2 

V>17 
& (Amax > 1.75 or 
V>19orV-R> 0.4) 

pfwsr > 0.6 
& rbcrossout < 0.75 

10 pts. > 2a 
& iVhi/A^pk > 0.9 

Xml-out/^dof < 4 
& Xrobust-out/^dof < 1.5 

None 

AxV(xLi/iVdof) > 300 
Amax > max(1.34, 1 + Aa) 
pfrdev2 < 90 



AXsN-ML > 
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Table 3 

Candidate Microlensing Events 



Event 


ID 


RA (2000) 


Dec (2000) 


V 


V-R 


^max 


t 


A 

^max 




la 


2.5628.5917 


05 


14 


44 


-68 


48 


01 


19, 


.75 


0. 


.56 


433.6 


34.4 


7.15 


1 


.083 


lb 


79.5628.1547 


05 


14 


44 


-68 


48 


00 


19, 


.73 


0. 


.55 


433.8 


34.2 


7.62 


1, 


.051 


4 


13.5961.1386 


05 


17 


14 


-70 


46 


58 


20, 


.25 


0. 


.14 


1023.4 


45.4 


2.92 


1, 


.415 


5 


6.5845.1091 


05 


16 


41 


-70 


29 


18 


21 


.15 


0, 


.76 


400.4 


75.6 


47.28 


1 


.512 


6 


7.7420.2571 


05 


26 


13 


-70 


21 


14 


19, 


.97 


0, 


.12 


573.6 


91.6 


2.43 


0, 


.763 


7a 


10.3802.872 


05 


04 


03 


-69 


33 


18 


20, 


.87 


0, 


.32 


840.0 


102.9 


5.91 


1 


.398 


*7b 


10.3802.494 


05 


04 


04 


-69 


33 


18 


19, 


.85 


0, 


.18 


832.9 


43.9 


1.63 


1 


.046 


8 


77.7307.4800 


05 


25 


09 


-69 


47 


53 


20, 


.31 


0, 


.21 


764.3 


66.4 


2.19 


1 


.733 




80.6468.2746 


05 


20 


20 


-69 


15 


11 


19 


.58 


0, 


.27 


976.4 


179.2 


1.95 


6 


.579 


lOa'^ 


1.3324.122 


05 


01 


15 


-69 


07 


33 


19, 


.60 


0, 


.23 


582.2 


41.8 


2.45 


1, 


.455 


lOb'^ 


18.3324.1765 


05 


01 


16 


-69 


07 


33 


19, 


.47 


0, 


.14 


581.8 


42.7 


2.38 


1, 


.455 


*ir 


11.8746.130 


05 


34 


21 


-70 


41 


07 


20, 


.04 


0, 


.85 


368.5 


280.7 


11.23 


2, 


.365 


12a'= 


11.8622.1257 


05 


33 


51 


-70 


50 


57 


21 


.34 


0, 


.25 


366.6 


137.4 


7.32 


1, 


.238 


12b'= 


14.8622.4762 


05 


33 


51 


-70 


50 


59 


21, 


.51 


0, 


.39 


365.7 


170.5 


7.31 


1, 


.237 


1 


80 7080 5384 


05 


24 


03 


-68 


49 


12 


21, 


.02 


0, 


.26 




1 00 1 


2 36 


1, 


.170 


14 


11.8871.2108 


05 


34 


44 


-70 


25 


07 


19, 


.37 


0, 


.01 


1767.8 


100.1 


3.37 


0, 


.751 


15 


10.4162.3555 


05 


05 


46 


-69 


43 


51 


21, 


.03 


0, 


.14 


1848.9 


36.8 


2.83 


0, 


.905 




79.4655.4621 


05 


09 


16 


-69 


08 


15 


18, 


.71 


0, 


.56 


1934.1 


27.0 


1.76 


0, 


.984 




9.5362.408 


05 


13 


35 


-70 


24 


43 


20 


.05 


0, 


.48 


2199.5 


24.1 


1.85 


1, 


.059 


18 


15.10554.465 


05 


45 


21 


-71 


09 


11 


19, 


.55 


0, 


.35 


1159.3 


74.2 


1.54 


1, 


.217 


*igc 


15.10669.178 


05 


46 


18 


-71 


31 


48 


19, 


,44 


0, 


.82 


687.5 


16.3 


1.41 


0, 


.734 


*20 


17.2221.1574 


04 


54 


19 


-70 


02 


15 


21, 


.35 


0, 


.57 


2151.2 


72.7 


2.95 


1, 


.298 


21 


17.2714.1058 


04 


57 


14 


-69 


27 


48 


19, 


.37 


0, 


.03 


589.7 


93.2 


5.64 


1, 


.620 


22 


22.5472.1126 


05 


14 


32 


-71 


09 


12 


20, 


.64 


0. 


.25 


1333.2 


229.3 


2.70 


1, 


.277 


23 


23.4143.256 


05 


06 


17 


-70 


58 


47 


20 


.27 


0, 


.31 


1138.8 


85.2 


2.41 


1, 


.492 


*24c 


24.2862.1187 


04 


57 


46 


-67 


41 


08 


20, 


.63 


0, 


.49 


2201.4 


186.5 


4.09 


3 


.927 


25 


24.3583.2286 


05 


02 


16 


-68 


00 


52 


19, 


.04 


0, 


.48 


1110.8 


85.2 


1.50 


0, 


.733 


*26'= 


47.1764.464 


04 


51 


11 


-68 


16 


41 


19, 


.46 


0. 


.21 


2170.6 


44.8 


1.88 


1, 


.305 


*27 


1.4289.1748 


05 


06 


35 


-69 


20 


48 


19, 


,24 


0, 


.07 


890.6 


50.5 


1.45 


1, 


.239 



The magnitudes and colors are fit baselines using the best available calibrations for each field as described 
in § H Time of peak magnification ^max is in JD - 2,448,623.5 (January 2, 1992). 

'^Events 1-12 appeared in A97. We number the current sample 1, 4 ... 27 to avoid any ambiguity with 
the previously published events. 

''Event 9 is the binary microlensing event; the parameters here are those resulting from a single-lens fit, 
and are not strictly appropriate. 

■^Probable supernova. 

*Events marked with "*" do not pass selection criteria A. 
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Table 4 

microlensing fits with blending 



Event 


^max 


t 


^max 


fov 


foR 


fv 


fR 





la 


433.6 


34.5 


7.19 


58.60 


93.30 





984 


1 


000 


1 


083 


lb 


433.7 


34.7 


7.83 


47.23 


77.38 





972 





982 


1 


052 


4 


1023.0 


83.3 


6.98 


40.85 


35.00 





322 





365 


1 


380 


5 


400.4 


109.8 


1.2e7 


15.98 


34.35 


1 


000 





457 





843 


6 


573.6 


92.0 


2.45 


42.63 


41.33 





981 


1 


000 





764 


7a 


840.1 


112.6 


6.87 


23.12 


25.51 


1 


000 





748 


1 


328 


*7b 


833.2 


46.5 


1.71 


61.37 


56.25 


1 


000 





761 


1 


042 


8 


764.3 


66.4 


2.19 


38.45 


32.52 


1 


000 


1 


000 


1 


735 




979.5 


143.4 




79.37 


83.30 





260 





170 


1 


755 


lOa*" 


582.1 


43.6 


2.56 


74.38 


71.49 


1 


000 





871 


1 


451 


1 Clh^ 




izo.o 


1 n 1 9 






n 
u 


lUO 


n 
u 


1 fin 


1 

1 


4:4: 1 


*11^ 


367.7 


436.9 


20.96 


31.24 


93.80 





998 





366 


1 


965 


12a^ 


367.1 


213.5 


13.41 


13.09 


15.59 





604 





440 


1 


184 


12b'' 


367.0 


1002.0 


65.94 


9.15 


12.45 





140 





091 


1 


138 


13 


1510.0 


222.7 


6.95 


17.79 


17.64 





219 





260 


1 


158 


14 


1768.0 


106.5 


3.67 


81.16 


59.87 





901 





874 





750 


15 


1849.0 


41.9 


3.48 


20.11 


16.88 





776 





735 





906 




1934.0 


27.0 


1.76 


148.00 


226.70 


1 


000 


1 


000 





985 




2200.0 


24.3 


1.88 


42.44 


61.96 


1 


000 





960 


1 


060 


18 


1159.0 


75.8 


1.58 


80.26 


91.15 


1 


000 





892 


1 


217 


*19b 


687.5 


18.6 


1.62 


86.28 


194.90 


1 


000 





465 





702 


*20 


2152.0 


99.4 


5.42 


12.30 


20.21 





628 





407 


1 


253 


21 


589.4 


141.5 


11.59 


91.01 


67.39 





499 





502 


1 


592 


22 


1333.0 


233.9 


2.73 


26.04 


26.31 





927 


1 


000 


1 


278 


23 


1139.0 


88.9 


2.61 


39.25 


41.98 


1 


000 





801 


1 


452 


*24b 


2201.0 


186.3 


4.09 


22.58 


32.85 


1 


000 


1 


000 


3 


941 


25 


1111.0 


85.3 


1.51 


112.50 


156.40 





924 


1 


000 





734 




2169.0 


260.0 


21.34 


82.81 


77.29 





056 





057 


1 


277 


*27 


895.0 


3247.0 


258.90 


85.17 


67.56 





002 





003 


1 


218 



Time of peak magnification t^ax is in JD - 2,448,623.5 (January 2, 1992). 

Tor the binary microlensing event (9) the fit parameters are given for the binary 
lens fit (Alcock et al. 1999b). Not all of the these parameters are appropriate for 
this fit. 

^Probable supernova. 

*Events marked with "*" do not pass selection criteria A. 
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Table 5 

SrpERNO\'A Type Ia Fits 



Event 


imax 


Uv 




6 


J yjv 


J UiX 


A / ' CLOJ 


la 


431.1 


37.44 


36.82 


0.75 


57.86 


92.39 


3.946 


lb 


431.1 


37.68 


37.10 


0.75 


46.74 


76.76 


3.283 


4 


1020.0 


38.11 


37.98 


0.75 


40.75 


34.88 


1.617 


5 


397.1 


36.94 


36.76 


0.75 


15.21 


34.57 


4.187 


6 


566.4 


39.49 


38.89 


-0.50 


42.99 


41.66 


0.865 


7a 


837.1 


37.44 


37.50 


0.70 


23.14 


25.60 


1.385 


*7b 


830.1 


38.50 


38.69 


0.75 


61.39 


56.26 


1.058 


8 


758.1 


39.63 


39.19 


-0.22 


38.51 


32.24 


1.750 


*9 


952.1 


39.02 


38.10 


-0.50 


84.93 


77.97 


7.030 


lOa*^ 


577.2 


37.72 


37.68 


0.59 


74.15 


71.18 


1.286 


lOb^ 


576.3 


38.02 


37.95 


0.33 


80.06 


67.04 


1.252 




362.5 


36.97 


36.26 


-0.42 


36.47 


98.99 


1.871 


12a^ 


363.2 


38.45 


38.24 


0.14 


13.33 


15.74 


1.101 


12b'' 


363.9 


38.66 


38.45 


0.11 


9.59 


12.77 


1.128 


13 


1497.0 


40.52 


39.56 


-0.50 


18.41 


18.15 


1.287 


14 


1758.0 


38.31 


38.06 


-0.50 


82.35 


60.47 


1.754 


15 


1845.0 


39.04 


39.10 


0.75 


20.30 


16.97 


0.948 




1930.0 


37.78 


36.91 


0.75 


148.00 


226.60 


0.976 




2196.0 


39.13 


38.71 


0.75 


42.26 


61.82 


1.036 


18 


1147.0 


39.49 


38.92 


-0.50 


80.91 


91.76 


1.406 




684.7 


38.52 


38.30 


0.75 


86.21 


194.80 


0.699 


*20 


2150.0 


38.89 


38.71 


0.72 


12.42 


20.49 


1.305 


21 


586.1 


36.10 


36.25 


0.75 


91.31 


67.57 


2.287 


22 


1318.0 


39.40 


38.89 


-0.50 


26.83 


27.32 


1.797 


23 


1128.0 


39.76 


39.22 


-0.50 


39.98 


42.12 


1.855 


*24a 


2183.0 


38.86 


37.96 


-0.50 


23.56 


34.24 


3.544 


25 


1103.0 


39.49 


38.49 


-0.50 


113.00 


157.80 


1.243 


*26* 


2165.0 


38.16 


38.06 


0.41 


83.10 


77.45 


1.222 


*27 


892.9 


:-!<S.(i4 


.S8.(i2 


0.7.1 


86.24 


GS.44 


l:.MY2 



Time of peak magnification ^max is in JD - 2,448,623.5. Events marked 
with "*" do not pass selection criteria A. 

''Probable supernova. 
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Table 6 

Supernova and Event Summary Table: 



Event 


XML/^dof 




^XSN-ML 


Galaxy? 


FoUow-Up 


Notes 


la 


1.083 


3.946 


6444 


No (HST) 






lb 


1.052 


3.283 


3685 


No (HST) 






4 


1.380 


1.617 


306 


No (HST) 


CTIO-l-spec. 




5 


0.843 


4.187 


4021 


No (HST) 






6 


0.764 


0.865 


173 


No (CTIO) 






7a 


1.328 


1.385 


71.0 


No (HST) 






*7b 


1.043 


1.058 


20.0 


No (HST) 






8 


1.735 


1.751 


28.0 


No (HST) 






*9 


6.269 


7.029 


1047 


No (HST) 




Binary ML 


10a'' 


1.451 


1.286 


-163 


Yes (HST) 






lOb** 


1.447 


1.252 


-226 


Yes (HST) 








1.966 


1.870 


-117 


Yes (CTIO) 






12a^ 


1.184 


1.101 


-113 


Yes fHST^I 






12b'* 


1.138 


1.128 


-12.0 


Yes (HST) 






13 


1.158 


1.287 


276 


No (CTIO) 


CTIO 




14 


0.750 


1.755 


1724 


No (HST) 


CTIO 




15 


0.906 


0.948 


45.6 


No (CTIO) 


CTIO 






0.985 


0.976 


-13.0 


Yes (CTIO) 


CTIO 






1.060 


1.036 


-31.0 


Yes (CTIO) 


CTIO 




18 


1.217 


1.406 


274 


No (CTIO) 






*;L9a 


0.701 


0.699 


-3.1 


Yes (CTIO) 






*20 


1.253 


1.305 


23.9 


No (CTIO) 






21 


1.592 


2.287 


428 


No (CTIO) 






22 


1.278 


1.797 


352 


? (CTIO 4m) 




ML parallax or rare SN? 


23 


1.452 


1.855 


263 


No (MACHO) 








3.941 


3.544 


-218 


Yes (MACHO) 






25 


0.734 


1.243 


232 


No (CTIO) 








1.278 


1.222 


-64.0 


No (CTIO) 




SN Type la or exotic ML? 


*27 


1.218 


1.302 


160 


No (CTIO) 







''Blended microlensing fits. 
^Probably supernova. 

This table summarizes the available information for each event, and a subjective 'ycs-or-no' determination of 
whether a background galaxy is present, along with the source of the image used in making the determination. 
Events marked with "*" do not pass selection criteria A. 
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Table 7 

microlensing events used, efficiency corrected t, and single event optical depths 











ri/10-9(A) 




Event 


t 






Tl/10-9(B) 



1 


34.2 


41.9 


44.5 


5.0 


3.8 


4 


45.4 


55.5 


59.0 


5.9 


4.5 


5 


75.6 


92.4 


98.1 


8.3 


6.7 


6 


91.6 


112.0 


118.9 


9.7 


7.9 


7 


102.9 


125.8 


133.6 


10.7 


8.7 


8 


66.4 


81.1 


86.2 


7.5 


6.1 


*9 


143.4 




143.4 




9.3 


13 


100.1 


122.4 


130.0 


10.5 


8.5 


14 


100.1 


122.4 


130.0 


10.5 


8.5 


15 


36.8 


45.0 


47.7 


5.2 


4.0 


18 


74.2 


90.7 


96.4 


8.2 


6.6 


*20 


72.7 




94.3 




6.5 


21 


93.2 


113.9 


121.0 


9.9 


8.0 


*22 


229.3 




297.8 




20.0 


23 


85.2 


104.2 


110.7 


9.1 


7.4 


25 


85.2 


104.2 


110.7 


9.1 


7.4 


*27 


50.5 




65.6 




4.9 



The quantity tst is the average actual event timescale for events in our Monte Carlo calcu- 
lations which are detected with an unblended fit timescale of t . For the binary event 9, the 
blended binary fit t value is used. The quantity ri is the contribution of each event to the total 
microlensing optical depth, computed using equation (|^). Columns are marked (A) or (B) to 
indicate which selection criteria was used in the efficiencies. 
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Table 8 

Total Model Independent Optical Depths 



Criteria & t # of events measured r(lO^) 9{degrees) 



A?-stat 13 l-lOtol 1.94 ± 0.29 

Af-fit 13 l-14lo:3 

A no <: -correction 13 0.94 

B ?-stat 17 l-29to:3 1-86 ± 0.23 

Bf-fit 17 l-24to:3 

B no f -correction 17 1.08 



The table entries show the microlensing optical depth r in units of 10"'', 
for the two selection criteria A & B using different t corrections (due to 
blending). The quoted errors are la standard errors computed as described 
in A97. The statistical t correction is preferred because it is unbiased. The 
'no t correction' values are for comparison only. 
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Table 9 

Optical Depth Confidence Intervals 



Selection Criteria 


# of events 






r(10-^ 


) for confidence level: 










0.025 


0.05 


0.16 


measured 


0.84 


0.95 


0.975 


A 


13 


0.60 


0.67 


0.83 


1.10 


1.47 


1.73 


1.86 


B 


17 


0.73 


0.81 


0.99 


1.29 


1.69 


1.97 


2.10 


A B average 




0.67 


0.74 


0.91 


1.20 


1.58 


1.85 


1.98 


A likelihood 


13 


0.40 


0.46 


0.62 


0.92 


1.32 


1.61 


1.76 


B likelihood 


17 


0.52 


0.58 


0.74 


1.05 


1.44 


1.73 


1.87 



The table entries show limits at various confidence levels on the microlensing optical depth r in 
units of lO"*" for different choices of selection criteria and different calculational methods. Rows 
marked A and B are model independent values found using the Monte Carlo method described 
in A97. The boldface row is the average of the model independent A and B calculations. The 
rows marked "likelihood" are background subtracted (that is, for halo microlensing only) and 
depend upon the model of the Galax y and LMC used, in this case model "S" with a dark LMC 
halo included, as described in § 3.2.3. 



Table 10 
Optical Depth Error Budget 



Cause 


Size (r/10-^ ) 


Relative Size 


Poisson 


0.4 


30% 


Exposure/Normalization of S/0 


0.25 


20% 


Selection Criteria 


0.25 


20% 


t bias 


0.05 


3% 


Binary star sources 


? 


? 


Exotic Lensing 


7 


7 
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Table 11 
microlensing by stars 



Population 


t(10-8) 


(t) (days) 


(O(kpc) 


r(io-V"^) 


^cxp (A) 


iVcxp (B) 


Thin disk 


0.36 


101 


1.3 


1.7 


0.38 


0.49 


Big thin disk (F) 


0.59 


101 


1.3 


2.7 


0.60 


0.79 


Thick disk 


0.20 


104 


3.6 


0.90 


0.20 


0.26 


Spheroid 


0.20 


129 


8.8 


0.90 


0.19 


0.25 


LMC disk(w/ halo) 


1.6 


120 


50 


5.8 


1.3 


1.7 


LMC disk(w/o halo).... 


2.6 


120 


50 


9.8 


2.2 


2.9 



This table shows microlensing quantities for various lens populations, with the density and velocity 
distributions and PDMF described in the text, r is the optical depth, (/) is the mean lens distance, 
and r is the total theoretical microlensing rate, but in all cases excluding bright lenses (see text). The 
expected number of events A^exp includes our detection efficiency averaged over the t distribution. The 
LMC values are averaged over the locations of our 30 fields. A^exp is the number of expected events using 
either selection criteria A (13 events), or criteria B (17 events). Two models of the LMC are considered, 
one with a dark halo and one without. Lensing from the LMC stellar disk only is shown in this table, 
lensing from the dark LMC halo is discussed elsewhere. 



34 



Table 12 
Maximum Likelihood Fits 



Events 


Coniment 


(Mo) 


/ML 


jML-iWH 

(101° Mq) 


(io-«) 


AT 

MW 


AT 

LMC halo 


AT 

stars 


S/13 


standard 




"■^^-0.07 


8.51^ 


101^3 


9.6 





3.0 


S/13 


standard 


54+0-26 








9.4 


1.1 


2.1 


S/17 


standard 




^•^^-0.08 




lll^ 


12.7 





3.9 


S/17 


standard 


r, 79+0.30 
' ^-0.20 


n OO+0.08 

"■^^-0.07 


9.111 




12.4 


1.4 


2.7 


B/13 


big halo 


0.68t°:i 


n i2+° °6 






9.7 





3.0 


B/13 


big halo 


U.OD_Q 22 


i2+° °5 

"■^^-0.04 




ntl 


9.8 


0.62 


2.1 


B/17 


big halo 


U-y^-0.28 


n u+° °^ 

^■^^-0.04 


lOtt 


ntl 


12.5 





3.9 


B/17 


big halo 


0.87t°:i 


14+0.05 

"■^^-0.04 


lOtt 




12.9 


0.78 


2.7 


F/13 


small halo 


^•-■^"-0.05 


n 'sn+°-22 

"■3^-0.18 


lOtt 




9.5 





3.2 


F/13 


small halo 


^•-■^^-0.06 


39+° " 

'-'•'->='-0.13 


8.01^3 


11-3 


7.0 


3.3 


2.3 


F/17 


small halo 


n 99+0.09 
'-'•^^-0.06 


0.57t°- 




nit 


12.5 





4.2 


F/17 


small halo 


n 25+°-i" 






ntt 


9.2 


4.3 


3.0 



The first column shows the model as defined in A96 and A97, and the number of microlensing 
candidates used; either 13 from selection criteria A, or 17 from criteria B. Model S is given by eq. ^ 
and has a typical size halo, Model B has a halo as large as possible, and Model F has a halo as small 
as possible with a large thin disk. Columns 3 & 4 show the maximum likelihood MACHO mass and 
halo fraction from Section 6.2.3| . Columns 5 and 6 show the implied total mass of MACHOs within 
50kpc of the Galactic center, and the resulting halo optical depth. For models with dark LMC halos, 
the sum of LMC and Milky Way MACHO optical depth is shown. Column 7 shows the number of 
expected events from the Milky Way Halo, column 8 shows the number of expected events from the 
LMC halo, and column 9 shows the expected number of events from stars (from Table 11). Every 
Milky Way model is shown twice, once with a dark LMC halo, and once with the dark LMC halo set 
to zero. See the text for more explanation. 
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Fig. 1. — An R-band image of the LMC, 8.2 degrees on a side (G. Bothun, private communication), showing 
the locations of the 30 MACHO fields used here. Also shown are the locations of the 17 microlensing candidates 
discussed in the text. 




Fig. 2. — Illustration of the cuts used to select microlensing candidates for criteria B. The x-axis is Ax^/ (x^j/A'dof ), 
where Ax^ = Xconst ~ Xmi improvement in x^ between a constant brightness fit and a microlensing fit. The 

y-axis is the fitted maximum magnification. The 29 lightcurves are shown as solid dots and arc labeled. The 
remaining symbols are explained in the figure and in detail in § 3.2. The solid lines show the final cuts for criteria 
B. The dotted line shows the same for criteria A. See text for details. 
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Fig. 3. — Illustration of the cuts used to select microlensing candidates for criteria B in the color-magnitude 
diagram. The 29 lightcurves are shown as solid dots and are labeled. The 'bumper' cut is outlined with solid lines 
and labeled with 'Amax > 1.75' for criteria B. The dotted lines shows the same for criteria A. The symbols are 
explained in the figure and in detail in § 3.3.1. The magnitudes displayed here use the rough global calibrations 
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Fig. 4. — The lightcurves for the 29 candidates (25 stars) in § B.2. 
show red and blue passbands. Flux is in linear units with la estimated errors 

brightness. Full lightcurves are shown with 2 day bi ni^ffig, insets of the event region s are unbinned. 



For each object, the upper and lower panels 
normalized to the fitted unlensed 

The full 



lightcurves can be found on the World Wide Web at [http : //wwwmacho . einu . edu . au/| The thick line is the fit to 
unblended microlensing (Table ^) , except for probable supernovae where both the blended fit (solid line) and type 
la fit (dashed line) are shown. 




Figure |^ (continued) 
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Figure |^ (continued) 
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Figure |^ (continued) 
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Figure |^ (continued) 
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Figure |^ (continued) 
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Figure |^ (continued) 
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Fig. 5. — Microlensing detection efficiency (normal- 
ized to Umin < 1) for the year 5.7 MACHO data, as 
a function of event timescale t . The solid line shows 
the 'photometric' efficiency computed for cut A, and 
the dotted line for cut B as described in § ^ For com- 
parison the corresponding curves from year-1 (A96) 
and year-2 (A97) are also shown. 



Selecfion Criferia A 




0.2 0.4 0.6 0.8 



Fig. 6. — The upper (lower) panel shows the cumula- 
tive distribution in Umin for set A (set B; excluding the 
binary event 9). The expected distribution is shown 
as a dashed line: a uniform distribution modified by 
our efficiency. The results of K-S tests comparing the 
observed and expected distributions are also shown. 
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Fig. 7. — Candidate microlensing events. The open circles indicate the position in the CMD of the blended MACHO 
object. The lines extend to the peak brightness of the event. The points correspond to the 200 MACHO objects 
nearest each event. The best available calibrations for each field described in § | have been used for this figure. 
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Fig. 8. — The upper (lower) panel shows the cumula- 
tive spatial distribution on the sky as measured from 
the optical center of the bar {a = 5h2Am, (5 = — 69°48') 
for set A (set B). The predictions of Gyuk et al. 
(1999) for the case of stellar LMC disk+bar self- 
lensing only {dotted line) and LMC disk+bar -|-halo 
self-lensing ( dashed line) are also shown. The mod- 
els of Gyuk et al. (1999) have been folded into our 
efficiency per field for each of the 30 fields reported on 
here. The results of K-S tests comparing the observed 
and expected distributions are also shown. 
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Fig. 9. — The contribution to the optical depth of 
eq. 1^ from events binned in timescale for selection 
criteria A (upper panel) and selection criteria B (lower 
panel). The solid line shows the observed values from 
the samples, with the number of events shown in each 
bin. The dotted curve shows the contribution to r 
which would arise from a single observed event with 
timescale t . The dashed line shows the total optical 
depth. 
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Fig. 10. — Measured microlensing optical depth is 
plotted as a function of the Wmin cut for selection cri- 
teria A [thick line) and selection criteria B [thin line). 
Our selected cuts for criteria A & B are marked with 
large dots plus 1 a error bars. The optical depth re- 
ported in A97 is also shown. 
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Fig. 12. — Likelihood contours for MACHO mass m 
and halo fraction / for Model S which has a typical 
size halo. See A96 for details of the model. The 
plus sign shows the maximum likelihood estimate and 
the contours enclose regions of 68%, 90%, 95%, and 
99% probability. The panels are labeled according to 
which set of selection criteria (A or B) is used, and 
whether or not an LMC halo with MACHO fraction 
/ is included. 
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Fig. 11. — Measured microlensing optical depth is 
plotted as a function of the Ax^ /{Xmif-^doi) cut for 
selection criteria A (thick line) and selection criteria 
B {thin line). Our selected cuts for criteria A & B 
are marked with large dots plus 1 a error bars. The 
optical depth reported in A97 is also shown. 
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Fig. 13. — Likelihood contours for MACHO mass 
m and halo fraction / for Model B which has a very 
large dark halo. See A96 for details of the model. The 
plus sign shows the maximum likelihood estimate and 
the contours enclose regions of 68%, 90%, 95%, and 
99% probability. The panels are labeled according to 
which set of selection criteria (A or B) is used, and 
whether or not an LMC halo with MACHO fraction 
/ is included. 



Fig. 14. — Likelihood contours for MACHO mass 
m and halo fraction / for Model F which has a very 
small dark halo, and a nearly maximal thin disk. See 
A96 for details of the model. The plus sign shows 
the maximum likelihood estimate and the contours 
enclose regions of 68%, 90%, 95%, and 99% probabil- 
ity. The panels are labeled according to which set of 
selection criteria (A or B) is used, and whether or not 
an LMC halo with MACHO fraction / is included. 
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